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Project N5ori-07856 
SURVEY 

In the advanced course which I have been giving this term, I have been treating 
the energy band theory. I had hoped to get further; but it has taken the entire term to 
cover the material. This includes just Chapter 1 of the proposed report on "Electronic 
Structure of Solids". Under the circumstances, it has seemed best to include just this 
chapter on energy bands in the Technical Report No. 4, and this report. Including also 
a bibliography of papers on energy bands, is being prepared for the printer, and should 
be sent out not long after this Progress Report. I had hoped also to treat the use of 
Wannier functions in discussing lattice imperfections and excitons, and various applica¬ 
tions of configuration interaction in solids, including ferromagnetism. Since this ma¬ 
terial was not covered, I shall propose to take it up in a continuation of the same course 
in the fall, and shall write the proposed additional chapters of the report on "Electronic 
Structure of Solids" at that time, bringing them out as a further Technical Report No. 5 
at a later date. 

This course on the theory of energy bands has furnished an occasion for exam¬ 
ining the relative merits of the various methods in current use for studying energy 
bands. The conclusion to which we have been led is that the orthogonalized plane wave 
method, perhaps modified to Herman's "augmented orthogonalized plane wave method", 
is the most promising method at present for accurate numerical work. Parmenter's 
attempt to apply the tight binding approximation accurately to chromium has shown the 
almost impossible complexity to which that method leads when we try to apply it accu¬ 
rately; the orthogonalized plane wave method seems clearly simpler. Parmenter ac¬ 
cordingly is looking into ways of applying this method to chromium. 

Tne great drawback of the orthogonalized plane wave method is that it is not 
readily adaptable to values of the propagation constant other than a few special points 
in the unit cell, such as the center of the Brillouin zone and a few other points. The 
reason is that at these symmetry points, a number of plane waves must have identical 
coefficients, so that the number of distinct coefficients is relatively small, even for a 
sum of a large number of plane waves, and we do not face a very difficult secular equa¬ 
tion. A method of interpolation between these symmetry points is highly desirable. 

With this in mind, it occurred to me that the tight binding approximation might furnish 
a good interpolation method, if the various integrals which occur in it were regarded as 
disposable parameters, rather than quantities to be calculated from wave functions and 
potentials. Dr. Koster and I have been examining this possibility, and it seems very 
encouraging. In particular, we have carried through a study of the diamond lattice, 
fitting the points determined by Herman by the orthogonalized plane wave method, and 
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our solution, which describes both the valence and conduction band, seems very likely 
to be rather close to the truth. We find that we can get this fit by using only nearest 
and second nearest neighbor interactions, and including no three-center or overlap in¬ 
tegrals; there are enough remaining integrals to use as disposable parameters. For¬ 
tunately, the integrals derived by fitting Herman's values have entirely reasonable 
values, having the correct signs and relative magnitudes, so that this procedure seems 
entirely reasonable. We are carrying this method further, considering somewhat more 
complicated crystals. This work Is report partly in the present Progress Report, part¬ 
ly in Technical Report No. 4. We are indebted to Dr. F. Herman, of R. C.A., for 
valuable help in this problem; he has also been independently considering the tight bind¬ 
ing approximation for interpolation purposes. 

The work on a two-electron model of ferromagnetism, on which Drs. Statz and 
Koster and I have reported in the two preceding Progress Reports, has been written up 
for submission to the Physical Review. Drs. Statz and Koster give in this Progress 
Report a two-dimensional model, coming a little closer to the truth than the one-dimen¬ 
sional model which I presented in the preceding Progress Report, and which is partly 
incorporated in the Physical Review paper. As for future work on this theory, we feel 
that we want to get closer to the actual structure of a ferromagnetic metal, such as 
nickel, and this demands knowing something about the Wannier functions and energy 
bands. We a»*e delaying work in this direction, on account of the feeling that the tight 
binding interpolation method, which I have just been describing, may well be adapted 
to such energy bands, and to the Wannier functions to be formed from them. Dr. Ki- 
kuchi, meanwhile, is looking at a very simplified two-atom model, in the hope that he 
will find a ferromagnetic ground state which may well have analogies to the ground state 
in more complicated systems. 

Dr. Pratt continues to work on the problem of antiferromagnetism, which is 
one of the most intricate problems which we face. His two communications, and one of 
Dr. Koster on an extension of Hund's rule, all have a bearing on this problem, whose 
status is something of the following sort. We have examined a number of models; and 
the difficulty which arises is that there is a great tendency for a ferromagnetic, rather 
than an antiferromagnetic, state to lie lower. This is probably a result of the method 
of approximation, but It illustrates the extreme sensitivity of the problem, and the need 
of examining it very carefully before we really understand the validity of the various 
approximations. In particular, Koster shows that Hund's rule, stating that the ferro¬ 
magnetic state will lie lower, holds for any case, molecular as well as atomic, in 
which our configuration interaction problem is an ordinary spin degeneracy problem 
between orthogonal orbitals. Probably as a result of this, Pratt finds that the model 
which he has set up of a triatomic molecule (MnOMn) ++ has a ferromagnetic state as 
its ground state, though we had hoped that we should find an antiferromagnetic ground 
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state, and so an explanation of superexchange. We are not satisfied, however, that 
this result is correct, for very likely a much more extensive configuration interaction 
is needed to lead to a possible antiferromagnetic ground state. 

The writer suspects that the actual situation is the following. He has pointed 
out on several occasions that the effective potential energy in which an electron of + 
spin moves may be different from that for an electron of - spin: on account of exchange 
effects, the electron of + spin has a lower potential energy when it is in a distribution 
of charge of other + electrons, and a higher potential energy when surrounded by elec¬ 
trons of - spin. If we then have an antiferromagnetic configuration of (MnOMn) + , the 
spins on one of the Mn ions (which will certainly all be parallel to each other in the 
ground state, by Hund's rule for the Mn ion) will point in one direction, say +, while 
those on the other Mn ion will point in the other direction, say -. An electron of + 
spin, then, will find a lower potential well in one Mn ion than in the other. This will 
hold for the electrons of oxygen, as well as for those of the manganese. Thus the elec¬ 
trons of + spin in the oxygen will find themselves effectively in a strong field, pulling 
them toward the manganese ion with ♦ spin, while the electrons of - spin in the oxygen 
will find themselves in an equal field in the opposite direction. These electrons will be 
polarized by these fields, just as if the fields were ordinary electrostatic fields, and 
this will reduce the energy of the molecule. A corresponding polarization will not occur 
in the ferromagnetic case, for then the electrons of each spin, in the oxygen, will find 
themselves effectively in zero field. The writer suspects that it is the lowering energy 
on accoimt of this polarization which stabilizes the antiferromagnetic state. 

This picture of the role of the oxygen atoms supplements the remarks made 
by the writer several years ago, regarding the alternating potential in which an elec¬ 
tron finds itself, in an antiferromagnetic crystal. It was pointed out that this alternat¬ 
ing potential would split the energy bands apart, the unit cell in the crystal being twice 
as great as ordinarily considered, and that in some cases we could well have the lower 
band filled, the upper one empty, thus explaining the non-conducting properties of cer¬ 
tain antiferromagnetic crystals which otherwise would be thought to be conductors. 

But this argument by itself does not show why the antiferromagnetic crystal has a lower 
energy than the ferromagnetic arrangement. Dr. Pratt reported some time ago on hav¬ 
ing looked into a one-dimensional model of such a system, and he found the ferromag¬ 
netic state to have the lower energy. It now seems likely that the presence of the oxy¬ 
gen is required to stabilize the antiferromagnetic state, in the way described above, 
though the splitting of the energy bands would take place just as described earlier by 
the writer. 

The problem is then simply to find the correct way to describe this polarize 
tion of the oxygens. This could be done by a configuration interaction with excited 
states of the oxygen atom, but this is a notoriously poor way to describe polarization, 
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from the point of view of convergence. It is much better to modify the wave function, 
in the way which Mr. Allen is investigating, using a single determinant made of orbit¬ 
als which are distorted as they would be in the presence of a field. In our case, the 
oxygen orbitals corresponding to electrons of one spin would be distorted in one direc¬ 
tion. those corresponding to the other spin in the other direction, so that they would no 
longer be Identical to each other. All this is easier to describe in simple language, 
than it is to carry out in detail. For one thing, we are dealing with a singlet state of 
the molecule, and we cannot say in any proper way that the spin of one manganese ion 
is pointed in one direction, that of the other ion in the other direction; the true wave 
function is a combination of wave functions of the two manganese ions, with a variety 
of orientations. The process of building up the wave function out of determinants is 
extremely complicated, on account of the many electrons involved; it is this part of 
the problem that Dr. Pratt is handling by his operator method. Finally, if the sus¬ 
picions regarding polarization are correct, it may be necessary to use non-orthogonal 
orbitals, and the technique for doing this properly, without neglect of the overlap in¬ 
tegrals. which may be the essential feature in the problem, is very complicated. The 
problems which Dr. Pratt is facing are not so much problems of visualizing the quali¬ 
tative nature of the solution, as those of formulating the solution in mathematical lan¬ 
guage in which we are sure that we are not throwing away the essential terms by mis¬ 
take. A simplified model which he has been considering, since treating the cases 
described in this Progress Report, gives hope for thinking that we are on the right 
track in our thinking about the polarization of the oxygen. 

The various investigations into molecular structure are progressing. An im¬ 
portant result of our consideration of configuration interaction has been the conviction 
that there are very few cases simple enough to handle by a complete configuration in¬ 
teraction. If we use only a partial configuration interaction, it then is very Important 
to have good molecular orbitals to start with. There are very few cases in the litera¬ 
ture in which really good calculations of molecular orbitals, by a self-consistent 
method, have been carried out. Accordingly, Dr. Meckler has considered the possi¬ 
bility of mechanizing the problem of determining self-consistent LCAO molecular or¬ 
bitals, and has found a way to do it on the Whirlwind Digital Computer, which he de¬ 
scribed in the present Progress Report. If this procedure can be actually put into 
operation, it might mean an important advance in the problem of finding molecular 
orbitals. 

The mechanization of the calculation of the integrals needed for the water 
problem, by Koster and Schweinler, is proceeding. This calculation, Meckler's on 
ethane, Kaplan's on ammonia, Barrett's on fluorine, and Corbato's on HF, all tie in 
together in the calculations of integrals, and general procedure, and are all progress¬ 
ing, though slowly. Unfortunately the calculation of HF by Mr. Merrifield, reported 
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on in preceding Progress Reports, had to be discontinued at a point where it had not 
yielded ureful numerical results; a continuation of this progiam is being studied by 
Mr. Corbato. Various other somewhat disconnected pieces of work are under way, 
as will be evident from the Progress Report. 

There has been one change in the personnel of the group since the preceding 
Progress Report: Dr. Statz has left the group to join the staff of the Raytheon Manu¬ 
facturing Company. Fortunately, he will still be in the neighborhood, so that he will 
be able to join in the general scientific activities of the group. Various other changes 
in personnel are in prospect for the summer and the next academic year. 


J. C. Slater 


1. A SIMPLIFIED TIGHT BINDING METHOD FOR ENERGY BANDS 


Among all the methods which have been used for studying energy bands, the 
tight binding approximation, or Bloch method, has one outstanding advantage: it gives 
simple analytical formulas for the energy as a function of the propagation constant. It 
has compensating disadvantages, however: to carry it through really rigorously, tak¬ 
ing into account all three-center integrals, overlap integrals, and other complicating 
features, is practically impossible, as Dr. Parmenter's attempt to apply it to the en¬ 
ergy bands of chromium has shown. It has occurred to the writers, however, that the 
method could be extremely useful as an interpolation method, for getting the energy as 
a function of propagation constant in a case where accurate values of the energy were 
known at certain symmetry points in the reciprocal space. This method has been writ¬ 
ten up in fair detail in the Technical Report No. 4, which will appear shortly, and for 
that reason only a sketch of the work so far done will be given in this Progress Report. 

Of all methods which have been used up to the present for energy bands, the 
orthogonalized plane wave method seems to the writers to be probably the one best 
adapted to getting accurate values of the energies. The study which Herman^ ^ has 
made of diamond and germanium gives one the feeling that the results are very reliable; 
particularly the indication that the method is converging rapidly as more and more 
terms are added to the series is very encouraging. On the other hand, the orthogonal¬ 
ized plane wave method is not well adapted for calculations except at particular sym¬ 
metry points in the reciprocal space. The reason for this is very simple: at a sym¬ 
metry point, many plane waves must have identical coefficients, so that the number of 
independent coefficients is far less than the number of plane waves which are being 
superposed to get an approximate solution of Schr&dinger's equation. Thus, for in¬ 
stance, the greatest number of independent coefficients which Herman used in his study 
of diamond was 16, so that he had to solve a 16-by-16 secular equation; but the wave 
function in this case was made of 146 plane waves. For a propagation constant not 
having special symmetry, we should still need 146 plane waves to get comparable ac¬ 
curacy, but now there would be no special relations between their coefficients, and a 
146-by-146 secular equation would be required, which of course would be beyond our 
present capabilities. 

We are then faced with a situation where we have very good energy values for 
certain values of the propagation constant, without an equally good way to interpolate 
between them. For instance, in the cases of diamond and germanium mentioned above, 
we have values for the center of the central Brillouln zone, and for the boundaries of 
the zone along the 100 and 111 directions. It has now occurred to Herman, as well 
as to the writers, that it might be possible to use a very much simplified tight binding 
method, but to choose the parameters entering into this method, not by direct calcula¬ 
tion from atomic wave functions and potentials, but by regarding them as arbitrary 
parameters, to be chosen so as to fit the values which were accurately known. We 
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might then hope that the solutions of the tight binding problem would give acceptable 
Interpolation values for the energy. Since we are using the method only for interpola¬ 
tion, we can afford to make simplifications which would be entirely inadmissible if we 
were using it for an original calculation: wc can disregard interaction integrals between 
distant neighbors, we can neglect all integrals except two-center integrals, and all over¬ 
lap integrals. We still keep enough terms, however, to have enough disposable con¬ 
stants to fit the values which we assume known; and of course we must keep such terms 
that we have a good hope that the calculation will agree with the correct one to an ade¬ 
quate accuracy. 

To illustrate the method, we have made a fairly complete study of the diamond 
crystal, fitting our constants to Herman's values. We start with the two interpenetrat¬ 
ing face-centered cubic lattices which form the diamond crystal. We form Bloch waves 
from s, p x> Py, p z orbitals on one of the lattices, and similar waves from the same 
orbitals on the other lattice. We have, then, eight Bloch waves, and we set up a secu¬ 
lar equation between these eight function for a given value of the propagation constant. 

In this secular equation, we disregard all overlap integrals between atomic orbitals on 
different atoms. The matrix components of energy involve integrals of the form 
J$ a *(r) H - R) dv, where and ^ are two atomic orbitals, (for instance, an s 

and a p x orbital), R is one of the displacement vectors from one atom to another, and 
H is the one-electron Hamiltonian operator of the periodic potential problem. It is 
well known that this integral can be reduced to a sum of integrals of the form J$^*lr) 

V(r - R') ^(r - R) dv, where the potential H is expressed as a sum of terms V(r - R'), 
spherically symmetrical potentials about the atoms located at positions R', to an ap¬ 
proximation which is sufficient for the present purposes. The integral above is a 
three-center integral, provided R and R' are different from each other and from zero. 

We disregard all such terms, and consider only the two-center integrals where R' 
equals R or zero. Then our integrals are like those found in a diatomic molecule. We 
can simplify them by expanding the $'s in terms of functions quantized with respect to 
the axis R of the pair of atoms. Thus if we are dealing with p orbitals, we have p? 
and pir components with respect to the axis. We have non-vanishing integrals only if 
both and ^ have the same component of angular momentum about the axis; that is, 
if both are a, or w, or 6 functions. We have integrals, then, of types which can be 
symbolized by the azimuthal quantum numbers of both orbitals (denoted by s or p in 
this case), and by the component of angular momentum about the axis (o- or it in this 
case). Such integrals can be denoted in an obvious way as (ssv), (spa), (ppcr), (ppv), 
and depend of course on whether we are dealing with nearest neighbors, next nearest, 
and so on; this can be denoted by subscripts 1, 2, etc. It is these various integrals 
which we regard as the disposable parameters of the problem. 

We have set up formulas for all matrix components of energy, diagonal and 
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non-diagonad, between the various Bloch functions in terms of these integrals. We 
then have, in our case of diamond, an eight-by-eight secular equation between the 
Bloch functions. This equation can be easily solved analytically at the required sym¬ 
metry points, and we can get solutions involving nothing worse than a biquadratic equa¬ 
tion along the 100 and 111 directions. We cam now try to fit all of Herman's values at 
the symmetry points. We find that, using nearest and second nearest neighbors, we 
do not have quite enough disposable constants to fit all of these values. However, we 
cam adjust the constamts to fit them all with only small errors. In other words, there 
are relations between Hermam's vadues which almost, though not quite, make tnem 
compatible with a solution involving only second nearest neighbors. When we do this, 
the resulting constamts aire Interesting, in that the integrals for second nearest neigh¬ 
bors come out much smaller tham those for nearest neighbors, suggesting that we 
readly have a convergent process. The integrals for nearest neighbors have the signs 
amd general magnitudes which we should expect; for instamce, (ppw), the interaction 
between pw orbitals on the two atoms, comes out considerably smaller numerically tham 
(pp<r), the interaction of p<r's on the two, as we should expect since the p<r's overlap 
more. 

When we use these constants, we then can get solutions for other points of the 
momentum space; we have carried out solutions along the 100 amd 111 directions (ex¬ 
cept that we have so far not solved for the states which require a biquadratic equation). 
The energies resemble the interpolations which have been used by Herman in his papers, 
but differ in some respects; these solutions are given in Technicsd Report No. 4. The 
quaditative behavior of the solutions is as given by Herman. Also the solutions in the 
neighborhood of the degeneracies, as in the neighborhood of the three-fold degenerate 
states at k = 0, have the behavior which we must expect from symmetry, amd which 
has been described by Shockley. ^ 

The method which we have described for diamond will be recognized as essen- 
tiadly equivadent to that proposed by Morita. ^ His formulas for the matrix components 
of energy ame essentially like ours, except that he does not express his integrals in 
terms of our quantities such as (ssa), etc. He finds his integrals from atomic orbitads 
amd potentials. We have worked backward from Morita's published values of his in¬ 
tegrals, to see the values which he must have found for the various integrals which we 
have determined by fitting Hermam's values. In doing so, we have observed what ap¬ 
peals to have been am error of some sort in Morita's integrals. Specificadly, he finds 
a numerically greater value for his integral which we denote by (ppw) tham for (pp<r), 
contrary to our findings, amd contrary to what one would expect from straightforward 
calculation. This apparent error in Morita's calculation seems to expladn amother fea¬ 
ture of his results. He found that the lowest energy in the conduction bamd, for k = 0, 
came from am s-like rather tham from a p-iike state, contrary to Herman's result. We, 
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however, have been able to find values of the integrals which not only put the levels in 
the same or«.’er as found by Herman, but agree quite well quantitatively with Herman's 
results. We feel, therefore, that the doubts about the tight binding approximation ex¬ 
pressed by Herman in the reference quoted are not justified. 

A different sort of tight binding approach has been carried out by Hall, ^ 
though he does not like to admit that his method is really a tight binding calculation. 

He starts with what he calls equivalent orbitals. We can describe these in the following 
fashion. We start with tetrahedral directed orbitals on each of the atoms of the crystal, 
directed toward the neighboring atom. We then make a symmetric combination of the 
two directed orbitals at the two ends of each bond, and call this symmetric combination 
an equivalent orbital. We have four such equivalent orbitals per unit cell of the crystal, 
so that we can make four Bloch combinations of them. We then can set up matrix com¬ 
ponents of the energy between these four Bloch functions, and can solve a four-by-four 
secular equation between them. This secular equation will have four roots, which may 
be expected to describe the four levels in the valence hand of diamond. It is this secu¬ 
lar equation which Hall has considered. 

The results which he has found are not in good agreement with those of Her¬ 
man; in particular. Hall, like Kimball and Hund and Mrowka^^ with the cellular approxi¬ 
mation, finds a doubly degenerate state whose energy is independent of the propagation 
vector, whereas Herman finds that the energy of this state is a function of k, and it is 
split in certain d J rections; our tight binding calculation agrees entirely with Herman in 
these respects. One way of describing this shortcoming of Hall's treatment is that he 
takes into account only interactions of nearest neighbors; if he considered more distant 
neighbors, he could get a suitable behavior. In fact, Herman^ has found it possible 
to get a good fit with his values by making up equivalent orbitals like Hall's, and using 
sufficiently distant neighbor interactions. The interactions must be assumed significant 
out to considerably more than the second nearest neighbors, however, which are the 
only ones which we need use in our method of formulation, in which we get the conduc¬ 
tion band as well as the valence band from and eight-by-eight secular equation. 

One can easily convert Hall's calculation into one equivalent to ours, and it is 
instructive to see how this is to be done. Just as Hall has set up four Bloch functions 
from the symmetrical equivalent orbitals, we could maka antisymmetric equivalent or¬ 
bitals, and make four Bloch functions from them. The four-by-four secular equation 
resulting from these would give a description of the conduction band. But we should 
find that there were non-diagonal matrix components of energy between the four Bloch 
functions formed from the symmetric orbitals, and the four formed from the antisym¬ 
metric orbitals. These non-diagonal matrix components vanish for k = 0, but become 
larger as we go out in the Brillouin zone. If we now consider all eight of these Bloch 
functions, and the eight-by-eight secular equation between them, we can show that it is 
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equivalent to the problem as we have set it up. We can get satisfactory results in this 
case by considering only nearest and second nearest neighbors, while to get satisfactory 
results from Hall's procedure we need terms coming from more distant neighbors. The 
easiest way to remove the zero-width bands, then, is to consider Interaction between 
the valence and conduction band. The point is that except at k = 0, the wave functions 
of the valence band are not constructed entirely from the symmetric combinations of 
orbitals on the two lattices, and the conduction band entirely from the antisymmetric 
combination; closer examination of the wave functions shows that the phase difference 
between the waves on these two lattices changes from zero or 180° at k = 0 (for the 
symmetric and antisymmetric solutions) to 90° at the edges of the band, a situation 
which can be described by Hall's method only with difficulty. 

We are examining the tight binding approximation further, and hope to come 
back to the relation between it and the Wannier functions for overlapping bands. 
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2. ENERGY BANDS IN CHROMIUM 


The original plan of calculating electronic energy bands in chromium, reported 
in previous Progress Reports, involved expanding the valence-electronic wave functions 
as arbitrary linear combinations of orthogonalized plane waves plus 3d Bloch functions. 
Each 3d Bloch function was to be a so-called Bloch sum of 3d atomic orbitals associated 
with the isolated chromium atom. In calculating matrix elements of the crystal poten¬ 
tial between two 3d Bloch waves, there arise terms involving three-center atomic in¬ 
tegrals (the one-electron integral of a product of an atomic potential centered on one 
lattice site and two 3d atomic wave functions centered on two different lattice sites). 
Because of the difficulty in accurately calculating three-center integrals, the above 
scheme is not feasible as it stands. 

Let us suppose, how ever, that the 3d atomic orbitals are modified in such a 
way that they vanish at distances greater than one-half the nearest neighbor distance in 
chromium. All three-center integrals will now vanish, and there will be no difficulty 
in calculating matrix elements of the crystal potential. Such a method should be general¬ 
ly useful among the transition elements. It has been suggested by Herman^^ that the 
same method (which he calls the augmented orthogonalized plane wave method) should 
be more rapidly convergent than the orthogonalized plane wave method alone when deal¬ 
ing with valence crystals (such as diamond) or semiconductors (such as germanium). 

With regard to transition elements, the motivation of the AOPW method is to represent 
by orthogonalized plane waves that portion of the crystal wave function resulting from 
both 4s atomic orbitals and the overlapping portions of 3d atomic orbitals and to repre¬ 
sent by Bloch waves that portion of the crystal wave function resulting from the non¬ 
overlapping portions of the 3d atomic orbitals. 

The method to be used in modifying the 3d atomic orbitals is more or less ar¬ 
bitrary, of course, but the following method, which has been used on chromium, seems 
physically reasonable and probably can be used whenever the AOPW method is applica¬ 
ble. Let us construct an atomic orbital by numerically integrating Schrddinger's equa¬ 
tion with a Hamiltonian containing the crystal potential, the numerical integration start¬ 
ing at any nucleus of the crystal. Since we desire an orbital of definite angular momen¬ 
tum (I = 2), we must use the spherical average of the crystal potential in Schrodinger's 
equation. This spherical average can be approximated by the isolated-atom potential 
V(r) for r < R and by the constant V(R) for r > R. The radius R is chosen such that 
V(R) » <V(r)>, where <V(r)> is the average value of V(r), averaged over the unit 
cell of the crystal. It seems reasonable to take the orbital equal to the isolated-atom 
3d orbital for r < R; i. e., we use the energy of the isolated-atom 3d orbital as the 
energy used in numerically integrating SchrOdinger's equation. For r > R, the numeri¬ 
cally integrated orbital will deviate sharply from the isolated-atom wave function since 
the spherical average of the crystal potential is lower than the isolated-atom potential. 
This will result in a radial node in the orbital at some value of r which is approximately 
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equal to one-half of the nearest neighbor distance in the crystal. (For chromium R = 

1. 55 atomic units, the radial node in the numerically integrated orbited occurs at 2.80 
atomic units, while one-half of the nearest neighbor distance is 2. 36 atomic units.) By 
cutting off the numerically integrated orbital beyond the radial node, we have formed 
a modified 3d orbital suitable for our purposes. 

Reference 

1. F. Herman, Ph. D Thesis, Department of Physics, Columbia University, January, 
1953, p. 40. 


R. H. Parmenter 
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3. A TWO-ELECTRON EXAMPLE OF FERROMAGNETISM 


In the previous two Progress Reports, J. C. Slater, ^ H. Statz, ^ and G. F. 
Koster^ have discussed a two-electron example of ferromagnetism. In this Progress 
Report we apply the same methods to a simple two-dimensional case. This case fur¬ 
ther confirms the necessity of degenerate bands for ferromagnetism and also illustrates 
some of the techniques which would be useful in a more detailed discussion of ferro¬ 
magnetism. 

A two-electron function in a crystal may be expanded in terms of products of 
Wannier functions associated with the various energy bands of the crystal. 


+(r,.r 2 ) 


(3.1) 


Here i refers to the band from which the Wannier function a^ arises, r^ and r^ are 
the coordinates of the two electrons, and R n is one of the primitive translations of the 
crystal. We know, however, that the symmetry properties of a wave function depend 
only on the invariances of the Hamiltonian and not on the number of electrons the wave 
function describes. We expect therefore that all the discussion of the irreducible rep¬ 
resentations of the space groups by Seitz^ and others^ will be applicable here. In 
particular, we are able to define a jf sector which expresses the fact that our two- 
electron wave function goes into <; K P limes Itself under translation through 
This allows us to write our U $ m ) in a simpler form 


iif- 


it ♦ it 


V ff n- O * ' 


F(i*„ - O 


(3. 2) 


More can be said since we know that our Wannier functions may have certain symmetry 
properties. Let us see how this influences our F's. We restrict our attention to K * 0. 
(The following arguments could be easily carried out for other if values.) We must 
first see what effect some operation a of our point group has on the wave function since 
we know that the wave functions for K * 0 must form irreducible representations of the 
point group of the crystal. It can easily be seen that 


*2> = t m F(q ’ 1 - 3 p ) aa j^**2 - ^p + 3 ) 

Here aa^(r) is the Wannier function which results from applying the operation a to 
aj(r). a" 1 is the primitive translation arrived at by applying a” 1 to ft m . Since we 
know the transformation properties of our Wannier functions from Eq. (3. 3) we are 
able to get relations between our F(^ m )’s. This is most easily seen by a simple example. 
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(A TWO-ELECTRON EXAMPLE OF FERROMAGNETISM) 


Assume that our two-electron wave function is derived from the Wannier functions of a 
single non-degenerate band for which aa(r) = a(r) is true for all o. Let us look for 
the two-electron wave function that forms a completely symmetric representation of 
our point group (a^r^, r^) = tjKrj, r^); for all a). By applying these relations to Eq. 
(3. 3) we get 


I F <ff m )a ( ?i - if p >.<? 2 • S p * S m ) = I F(«- 1 S m ) a(r, - S p )a ( r 2 - R p ♦ Rj 


therefore F(it ) = F(a" 1 it ) 


(3.4) 


The additional symmetry property of our two-electron wave function concerns its value 
of total spin. This is determined by the relation 

F(-S m ) = tF(S m ) (3.5) 


The - sign refers to the triplet and the + sign refers to the singlet. 

Let us now apply the prece -ing discussion to a two-dimensional square lattice. 
We wish to determine our F's by minimizing the expectation value of a two-electron 
Hamiltonian 


II = Hj + H 2 + g 12 (3.6) 

where Hj and H 2 are the one-electron Hamiltonians of our two electrons in a periodic 
potential and gj 2 is the Coulomb repulsion between two electrons or two holes. We 
study first a non-d«»generate band for which aa(r) = a(r) for all a of the point group. 
This gives rise to the secular equation (Ref. 1, fiq. (1.6)) 


x 2*(it ) F(it n + it ) + 6(it n , 0) V Q = EF(it n ) ( 3 . 7 ) 

P H 


Here 

/(itp) - /atf;) Hj a(r t - K p ) 

Mr ) a(r.) a(r,) i(r 2 ) 

V o = 2J — 1 - - --- - • (aalaa) 

r i2 
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(A TWO-ELECTRON EXAMPLE OF FERROMAGNETISM) 



Fig. 3-1 

Square lattice in Rj - R 2 space. The figure gives the notation 
for the g lattice points and the various "planes" of reflection. 

In this secular equation the assumption has been made that two electrons interact only 
when they are at the same lattice site. We also assume ^(ff p ) a 0 except for nearest 
neighbor interactions. For boundary conditions in relative coordinate space we make 
the assumption that F(lt n ) vanishes on the sides of a square that contains nine lattice 
sites. Fig. 3-1 contains a diagram of this region of relative coordinate space as well 
as labels for the various lattice points under consideration. Table 3-1 gives the vari¬ 
ous irreducible representations of group C^ y (the point group of a square lattice). The 
a's are reflections through the "planes" indicated in Fig. 3-1. C 4 is a clockwise rota¬ 
tion through 90° and = C^ 2 . 

Table 3-1 

Character table for the group of the square lattice C^ y 
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(A TWO-ELECTRON EXAMPLE OF FERROMAGNETISM) 


We will now discuss in detail one of the symmetries of twc-electron wave 
functions for K = 0 under the above assumptions. Let us look at those wave functions 
with Aj symmetry. Using the relation (3. 4) we see at once that 

F( 1, 1) = F( 1, - 1) = F(- 1,-1) = F(- 1,1) = * 

F(l, 0) = F(0, 1) = F(- 1, 0) = F(0, - 1) = | (3. 8) 

F(0, 0) = a 


Substituting these relations into secular equation (3. 7) we find 


therefore 


a(V Q - E) + |(8/) + 1(0) = 0 
a(2<f) + |(-E) ♦ 1(4/) = 0 
a(0) + |(4/) + 1(-E) - 0 

E 3 - V E 2 - 32/ Z E + 16/ Z V = 0 

o o 


(3.9) 


(3. 10) 


Here / is the nearest neighbor interaction (that is, the matrix element of the one- 
electron Hamiltonian between a Wannier function at (0, 0) and one at (0, 1)). The equa¬ 
tion (3. 10) can be solved for small / giving the result 


E 




V 

o 


(3.11) 


These states are clearly singlet states since F(R r ) 3 + F(-^ n ). A similar discussion 
can be carried out for the other symmetries. First we write down the relations be¬ 
tween the F's which result from the irreducible representation which we desire and the 
symmetry properties of the Wannier functions. We then use these relations to simplify 
the secular equation (3. 7) and finally we find approximate eigenvalues of the secular 
equation if these equations are above .ie second order. 

For A^ symmetry we find that all the F*s must vanish and therefore we cannot 
find a state of this symmetrv. 
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(A TWO-ELECTRON EXAMPLE OF FERROMAGNETISM) 


For Bj symmetry we find 

F(l.l) = F(-l.l) = F( 1, - 1) = F{- 1,-1) = 0 
F(0, 1) = F(0, 1) = F( 1.0) = - F( -1, 0) = | 
F(0, 0) = 0 

which leads to the secular equation 


|(-E) = 0; E = 0 


This is a singlet. 

For symmetry we find 

F( 1, 1) = F(- 1, - 1) = -F(- 1. 1) = - F( 1, - 1) = | 
F(0, l) = F( 1,0) = F(-1.0) = F(0, - 1) = 0 
F(0, 0) = 0 


This gives rise to the secular equation 


|(-E) = 0; E = 0 


(3. 12) 


(3. 13) 


(3. 14) 


(3. 15) 


once again this is a singlet. 

For E symmetry we have a doubly degenerate wave function. In this case we 
set up two wave functions 

+« ' I**<*m> » (r i ' ff p> a( ~2 - S p * B m> 

*j * E F3 ' <B m> “ (? 1 • S p> al? 2 • B p * S m> 

Here 4» x and transform as x and y do under the operations of the group of the 
square. Once again this gives us relations between the F's 

F x (l, 1) = F x (- 1, 1) = -F x (- 1,-1) = - F x (l, - 1) = | 

2 (3.17) 

F^l, 1) = -F^- 1, 1) = -F y ( - 1. - 1) = F^l, - 1) * | 
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(A TWO-ELECTRON EXAMPLE OF FERROMAGNETISM) 


F X (0, I) = - F x (0, - I) = F*(1,0) = -F^-l.O) = -2- 

F X (1,0) «= F x (- I, 0) = F^(0, I) = 1^(0, - 1) = 0 (3. 17) con'd. 

F x (0,0) = 1^(0, 0) = 0 

The secular equation 'or the x-like functions is 

“ (- E) ♦ £(4,?) = 0 

2 (3. 18) 

-A-(2/) + £(- E) = 0 
V2 2 

with roots 

E = t ijis 

Thus we get two doubly degenerate eigenvalues (a total of four wave functions). These 
states are all triplets (F(^ R ) - - F(- i$ n ))* Fig. 3-2 contains a graph of energy as a 


*UI 



Energy as a function of nearest neighbor inter¬ 
action £ for the case of non-degenerate bands- 

function of £ for the two-electron states arising from a non-degenerate band. We no¬ 
tice that the singlet state of Aj symmetry is the lowest state. The fact that the singlet 
lies lower than the triplet we interpret as an absence of ferromagnetism. 

We now consider the case of two degenerate bands which we describe in terms 
of Wannier functions with x- and y-like symmetry. ^ We will call them a(r) and b(r) 
and their symmetry properties are given in Table 3-2. 
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(A TWO-ELECTRON EXAMPLE OF FERROMAGNETISM) 


Table 3-2 

Symmetry properties of Wannier functions a and b 
Operation b a 

E b a 


dl 

*d2 


- a 
b 

- b 


- a 

- b 


We write our two-electron function as 


♦*'l' -V ■ I F aa (S m> a(? l - B p> a(? 2' B p * B m> 

+ £ F bb< B m> b < ? l ' B p> b ( r 2 - B p * B m> 

* X F „b< B m> a < F l - B p> - B p * B m> 

* 3>ba (B m> b < b l - B p> a < b 2 ' B p * B m> 


(3. 19) 


If we assume that there is only interaction between electrons when they are at the same 
lattice site, these functions give rise to the secular equations^ 

I *'a (B p> F aa (B n * B p> * MB n'°> [ F .. (0 > (aalaa) + F b b (0 > (ab ' ab >] * EF a a (B n> 

X *'b (B p> F bb< B „ * B p> * !(5 n .0)[F bb (0)(bblbb) * F aa (0)(ablab)] = EF bb (B n ) 

S ['a (B p> * 'b (B p>] F ab' B „ tB p) * « B n’°>[ F ab‘ 0 >< aalbb > * F ba «»U b l» b >] = EF ab < B n> 

l ['a (B p> * * b (B p>] F ba (B n + B p> + « B a’ °> [ F ba<°> (aa,bb > * F ab‘° Bablab >] = EF ba (B n> 

P (3.20) 

The secular equations break into two parts, one connecting the F &a and F fab type two- 
electron wave functions and the other connecting the F ab and F^ a types. This splitting 
is due to our approximation in considering interactions of electrons only when they are 
on the same lattice site. In what follows we make use of these abbreviations 
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(A TWO-ELECTRON EXAMPLE OF FERROMAGNETISM) 


V Q = (aalaa) = (bblbb) = J a^) a^) g J2 a^) a(r 2 ) d-Tj dr 2 

J = (ablab) = /a^) b(r“ 2 ) g 1£ a^) b(r z ) d T l dT £ 

K = (aalbb) = /afr^) afr^) g J2 atr^) b(r 2 ) dr j dr £ 

We also make the assumptions that only those Fourier coefficients of the one-electron 

energy that represent nearest neighbor interactions are non zero. We call 6^ the in¬ 
teraction between an x-like Wannier function and one displaced a unit translation in the 
x direction. We call £ z the interaction between an x-like Wannier function and one 
displaced a unit displacement in the y direction. A new feature comes into this calcu¬ 
lation of the energies of the two-electron functions associated with the degenerate bands. 
In this case it sometimes occurs that the symmetry of the wave function is not sufficient 
to determine whether the wave function is a singlet or a triplet. When this is the case 
there can be both a singlet and a triplet of the same symmetry. Instead of presenting a 
detailed discussion of the many symmetries and multiplicites we present the results in 
Table 3-3. The arguments that go into deriving the tabulated results are exactly analo¬ 
gous to those already presented for the non-degenerate band. The first column (I) rep¬ 
resents the multiplicity of the state. The second column (II) gives the symmetry. The 
third column (III) tells whether the wave function comes from the F fta type secular 
equations or the F ab type secular equations (see Eq. (3. 20)). In the fourth column (IV), 
there is listed the secular deterrr nant and in the final column (V) is listed the approxi¬ 
mate eigenvalues for £ ^ and £ z small. In this column we have used the abbreviation 
that £^ = c £ z . 

For the case of a two-electron wave function drawn from two degenerate bands 
it is not possible to draw an unambiguous conclusion as to whether the triplet or the 
singlet lies lower as one can see from Table 3-3. Of those states in which one electron 
is in each band (F &b type), it is clear that the lowest state of A 2 symmetry is a triplet 
lower than any other state of the F &b type. Of the wave functions in which both electrons 
are in the same band (F &a type), the lowest state is a singlet of Bj symmetry. As to 
the competition between these two, the lowest one will be determined by the actual 
values of the parameters involved. Either one of these two could lie lower. Thus we 
see that in the case of non-degenerate bands the singlet lies lower whereas in the case 
of degenerate bands it is possible that the triplet may lie lower. 

References 
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4. A SYSTEM OF TWO p ELECTRONS 


In connection with the work on the two-electron example of ferromagnetism 
developed by other members of the Group, ^ there occurred a question regarding the 
energy levels of a system of two Ni atoms as a function of the internuclear distance; 
the reason is that the properties of the pair of atoms might be supposed to give infor¬ 
mation about the ferromagnetic substance in bulk. 

In order to obtain qualitati /e information about this system, particularly 
about the multiplicity of the lowest lying state, the problem was simplified several 
steps as follows. A Ni atom which has one hole is 3d orbits is first reduced to an atom 
with a hole in p orbits, because both d and p levels give rise to degenerate bands. 

A hole in p orbits behaves qualitatively similar to an electron in p orbits. Further¬ 
more, the Gaussian type of wave functions may be enough to give the qualitative infor¬ 
mation. 

After these considerations, now the energy levels of a system of two atoms 
each having one p electron is being calculated using Gaussian wave functions, as a 
function of internuclear distance, and taking into account the configuration interaction. 

Reference 

1. H. Statz, Quarterly Progress Report, Solid-State and Molecular Theory Group, 

M.I.T. . January 15, 1953, p. 23; G. F. Koster, ibid. p. 31; J. C. Slater. 

April 15. 1953, p. 6. 
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5. SPIN OPERATOR METHOD -- APPLICATION TO CONFIGURATION INTERACTION 


A detailed account of an operator method of constructing eigenfunctions of S 2 
was presented in the last Progress Report. The application of this technique to con¬ 
figuration interaction has proved rather successful and some of the details of its use 
will be given here. The following remarks will be largely confined to the treatment of 
the Mn— ——Mn ++ model of antiferromagnetism discussed elsewhere in the present 
Progress Report. 

It was shown in the April 15 Quarterly Progress Report that there exists an 
operator O when operating on a 2N electron determinant whose M g value is zero, 
produces a linear combination of 2N electron determinants such that the resulting state 
is a singlet. The only condition imposed was that all the one-electron spatial orbitals 
be distinct, but not necessarily orthonormal. It was also pointed out in that Report 
that this operator method could be used to set up the matrix components of the Hamil¬ 
tonian between singlets corresponding to different configurations. Suppose Pj is a 
product function of space-spin orbitals 

Pl = Ul (l) q( 1) u 2 (2) q(2) ... u n (N) a(N) u N + ^N ♦ 1) p(N ♦ 1) ... u 2N (2N) p(2N) (5. 1) 

of total M g = 0. Then one of the singlet states that can be formed from this collection 
of spatial orbitals is 


4l» = AOpj (5.2) 

where O is the spin operator and A is the antisymmetrizing operator. The matrix 
component of the energy between two configurations for this particular type of singlet 
is 


<+! H +2> = <AOpj H AOp 2 > (5.3) 

We now make use of the following properties of the spin operator O: (1)0 commutes 
with a spin-free Hamiltonian. (2) O commutes with the antisymmetrizing operator A. 
(3) The square of O is equal to a constant times O. These properties allow one to 
write Eq. (5. 3) as 

<4#!|HI4» 2 > = (N + 1) 1/2 <A (2N) 0 UN) P jlHlp 2 > (5.4) 

where pj and p 2 are spin product functions corresponding to different configurations 
and A^^) and signify that these operators are functions of the coordinates of all 

of the electrons. 

In the ground configuration of the Mn * - —Mn ++ system the orbitals con- 
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(SPIN OPERATOR METHOD -- APPLICATION TO CONFIGURATION INTERACTION) 


sidered are the five Mn 3d orbitals each singly occupied and the oxygen 2p orbitals 
which arc assumed to hold six electrons. Thus we have a sixteen-electron system but 
only thirteen distinct orbitals. The object of these remarks is to generalize the spin 
operator O in order to handle systems composed of spatial orbitals all of which are 
not distinct. By doing this we will be able to cast matrix elements in the convenient 
form given in Eq. (5. 4). 

In the investigation of the Mn^i-O^—Mn ++ system the singlet considered for 
the ground configuration was formed by combining two S = 5/2 states made up from the 
Mn 3d orbitals to give a resultant ten-electron singlet. The remaining six electrons 
were assigned to the three 2p oxygen orbitals. In terms of the spin operator this six¬ 
teen-electron singlet is written 

$ = A (l6) [«,(!) o( 1) Uj(2) 0(2) . . . u 3 (5) «(5) Uj (6) 0(6) • O (l0) {u 4 (7) «(7) u & (8) a(8) 

u 6 (9)a(9)... u g (ll) u(ll) u 9 (12) 0(12) u, 0 (13) 0(13) ... u |3 (16) 0(16)}] ^ ^ 

where u ^, u^, and u 3 represent the 2p functions and u^ through u g are the 3d functions 

on one Mn ++ ion and u Q through u. , are the 3il orbitals on the other Mn ++ ion. An es- 
' D (10) 

sential point appears here; that is the ten-electron spin operator O' ' does not com¬ 
mute with the sixteen-election antisymmetrizing operator A^ l6 \ Therefore, unless 
this singlet can be set up using the full sixteen-electron spin operator the ma¬ 

trix components cannot be haidled as in (5. 4). 

Consider a 2N electron system with 2(N - n) distinct orbitals which will be 
kept singly occupied and n distinct orbitals which will be doubly occupied. The or¬ 
bitals may be catalogued in the following manner 



Group 1^ 




Group 


Group Ig 


"No 


Group Ilg 


(5.6) 


If all of the functions were distinct, the 2N electron singlet under discussion would be 
written 

* = a (2N) 0 (2N) {u 1a ( 1) a(l) ... u N (N)«(N)u 1b (N+1)P(N+1)... u^UN) 0(2N)}(5. 7) 
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(SPIN OPERATOR METHOD — APPLICATION TO CONFIGURATION INTERACTION) 


Suppose now we set Uj^ equal to uj^ for j going from 1 to n in the singlet 
given in (S. 7). Such a change cannot affect the multiplicity of the total wave function. 
However, many of the determinants which comprise (5. 7) will vanish as a result of 
assigning the same spin to orbitals whose spatial parts are now assumed identical. Of 
the non-vanishing determinants, many will differ from others merely by interchanges 
of rows and hence these determinants can be combined. The net result of this change 
of spatial orbitals is that the normalization of the resultant singlet is destroyed. O' 21 ^ 
must be modified so that the normalization can be preserved under such a change. 

It can be shown that the necessary modification of O^ 2 *^ is to multiply it by a 
constant factor which depends only on the total number of electrons, 2N, and the num- 
be of identical orbitals, n. This factor is 


[ 


n! (N - k)«1 - >/ 2 
(n - k)! N!j 


(5.8) 


Once the matrix elements have been reduced to the form of Eq. (5. 4), O^ 2 *^ 
can be replaced by an equivalent operator for the purpose of taking the matrix com¬ 
ponent. The equivalent operator is obtained by omitting entirely the one-electron op¬ 
erators S. + , Sj” where i and j refer to electrons occupying an orbital which is doubly 
filled. More precisely, this operator is 0 2 ^ N " and it operates only on the spin co¬ 
ordinates of the electrons in groups 11^ and H Q . This is a valuable result as it ma¬ 
terially reduces the amount of work necessary to find the matrix components. There¬ 
fore, in the case of Mn—O—Mn the matrix component of the energy between the singlets 
corresponding to the I th and j th configurations can be expressed as 


<^lHl^>= (N - n + 1) 1 / 2 <A (16) O {10) p t IH| Pj > (5.9) 

where N is 8 and n is equal to the number of doubly filled orbitals which is three. 


G. W. Pratt, Jr. 
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6. ANTIFERROMAGNETISM 


At the present time attempts to give a quantum mechanical description of anti- 
ferromagnetism have been largely restricted to the treatment of simple models which 
supposedly would show a preference for spins to be oriented antiferromagnetically. A 
configuration interaction treatment of the Mn^i——Mn ++ system is reported here. 

The Mn—Mn distance as determined by x-ray scattering is 4. 426 A°^ and it seems 
unlikely that the Mn 3d orbitals located about centers this distance apart would overlap 
appreciably. This overlap was neglected in this investigation. By considering excited 
configurations in which there is an unpaired spin on the oxygen, a spin dependent coupl¬ 
ing will appear which causes the non-overlapping charge distributions associated with 
the Mn ions to interact through the oxygen. This is commonly referred to as super- 
exchange interaction. ^ The object of this work was to examine the nature of this 
coupling to determine whether it would tend to orient the spins of the Mn electrons 
ferromagnetically or non-magnetically. 

The ground configuration was taken to be two Mn ++ ions and an O - ion. The 
Mn * ion has five 3d electrons outside closed shells and the O - ion was taken to have a 
rare gas configuration. The excited configurations considered were those obtained by 
allowing an electron in an orbital associated with the O' ion to occupy one of the Mn 3d 
orbitals. The only symmetry considered in the work was of the £ type as the ground 
state has this symmetry. 

The one-eleciron functions used were assumed to be localized and orthonormal 
although the results obtained are independent of their exact nature. Table 6-1 shows 
the assignment of electrons for the various configurations. The Mn 3d orbitals are de¬ 
noted by: mj = 2 (6 + ), m j = 1 (w + ), m j = 0 (a), nij = -1 (*_), and = - 2 (6 _). 

The oxygen orbitals are m } = 1 (P + ), nij = 0 (P Q ), nij = - 1 (P_). 

The states of maximum multiplicity, S = 5, for each configuration and the 
corresponding secular equation was set up. The energy of the lowest S = 5 state was 
compared with the energy of the lowest singlet found by constructing one particular 
singlet for each configuration and examining the resulting secular equation. Let us con¬ 
sider the choice of this singlet state. 

For each configuration there are ten singly occupied orbitals and, therefore, 

42 orthogonal singlets. From neutron diffraction studies of MnO, it is thought that the 
Mn ++ ions can be divided into two sublattices each sublatHce being ferromagnetically 
coupled:^ the singlet state most closely resembling this situation is that formed by 
combining the two sublattices which are separately in a state of maximum spin to give a 
resultant singlet for the entire system. The analogue of this picture in the model being 
considered here is for the two Mn ions to be each in a state of maximum multiplicity 
and for these two S = 5/2 systems to combine to give a singlet ftr the ground configura¬ 
tion. The singlet for each excited configuration was formed by the combination of two 
S = 5/2 systems. As it is only possible to form one independent singlet in this manner. 
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Table 6-1 



for each configuration one S = 5 and one S = 0 state was considered. 

The states which are even with respect to a reflection in a plane perpendicular 
to the symmetry axis and passing through the oxygen are: 

§ 8 s |b 
*0 M 

•i 8 a *z ♦ +3 

»2* * +4 * +5 

- (»• 1 ) 
*3 S * +6 * +7 

«4 8 ’ + 8 * +9 

» 5 8 * +10 * +11 


From these six states the states which have the + symmetry with respect to reflection 
in a plane containing the symmetry axis of the model can be formed. These four states 
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*3 = # 2 g + # 3 g = +4 + +5 + +6 + *7 

♦4 = # 5 g = +10 + +11 


(6. 2) con'd. 


The S = 5 and S = 0 states corresponding to these four £ + states were set up 
and the corresponding secular equations written down. The S = 5 states were handled 
by conventional methods in that the states ^(S = 5) are just a single determinant. 

The singlet states »J^(S = 0) are, however, much more complicated as these singlets 
are each expressed as a linear combination of 252 determinants. It would be an enor¬ 
mous task to compute the matrix components of the energy for these states by the cus¬ 
tomary procedures. These matrix components were found by the spin operator method 
reported in the preceding Quarterly Progress Report. The essence of this method is 
that one can construct a state of definite multiplicity by operating on a single determi¬ 
nant of one-electron space-spin orbitals with an operator which is essentially a poly¬ 
nomial the elements of which are step-up and step-down spin operators taken to vari¬ 
ous powers. This operator has the properties that it is hermetian, it commutes with 
a spin free Hamiltonian, it commutes with the S Z operator, it commutes with the anti¬ 
symmetrizing operator, and its squa~e is equal to a constant times itself. These prop¬ 
erties allow one to use the operator to take matrix components of the Hamiltonain be¬ 
tween singlets corresponding to all configurations. ^ Some of the details involved are 
given in another section of the current Progress Report. 

In terms of the spin operator O and the antisymmetrizing operator A, we 
can write the four 1 Y states as 

♦ j(S = 0) = AOpj 

♦ 2 (S = 0) = ACXp 2 + p 3 + Pg + p 9 ) 

(6. 3) 

4> 3 (S = 0) = AO(p 4 + p 5 + p 6 + p 7 ) 

♦ 4 (S = 0) = ACXp l0 + p u ) 

where the p^ are spin product functions referring to the various configurations. 


P t = 6 + (l) a(l) 6 (2) a(2) * + (3) a(3) w_(4) a(4) <r(5) a(5) P q ( 6) 0(6 ) P q (7) 0(7) P + ( 8) a(8) P + (9) 0(9) 
P.(10) a(l0) P_( 11) 0(1 l)«r'(l2) 0(12) v_'(l3) 0( 13) w + '(14) 0( 14) 6_'( 15) 0( 15) 6 + '(16) 0(l6) 

(6.4) 
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p 2 = 6 + ( l) a ( i) . . . a{5) a(5) P q ( 6) a(6) P q ( 7) p(7) P + (8) q( 8) P + (9) p(9) P_( 10) a( 10) 6_( 11) p(l 1) 
<r'( 12) p( 12) . . . 6 + '(l6)p(l6) (6.5) 

p 3 = 6 + ( 1) a( 1) . . . <r(5) a(5) P q ( 6) a(6) P q ( 7) p(7) P + (8) a(8) P + (9) p(9) 6_'( 10) a( 10) P_( 11) P( 11) 
*'(12)P(12) ... 6 + '(l6)p(l6) (6.6) 

p 4 = 6 + ( 1) a(l) . .. a(5) a(5) P q ( 6) a(6) P q ( 7) p(7) P + (8) a(8) P + (9) p(9) P_(10) a(l0) *_(l 1) p(l 1) 
*'(12)P(12) ... * + '(l6)p(l6) (6.7) 

P 5 = 6 + (l)a(l) ... <K5)Q(5)P o (6)a(6)P o (7)p(7)P + (8)a(8)P + (9)p(9)tr_'(l0)a(l0)P_(ll)p(ll) 
o-'(l2)p(12) ... 6 + '(l6)p(l6) (6.8) 

p 6 = 6 + (l)a(l) ... «K5)a(5)P o (6)o(6)P o (7)p(7)P + (8)a(8)ir + (9)p(9)P_(l0)a(l0)P_(ll) p(ll) 
«r'(l2)p(l2) ... 6 + '(l6)p(l6) (6.9) 

P? = 6 + ( 1) a(l) . . . «r(5) a(5) P q ( 6) a(6) P q ( 7) p(7) w + '(8) a(8) P + (9) p(9) P ( 10)a(10)P (11) p( 11) 

«r'(l2)p(l2) ... 6 + '(l6)p(l6) (6. 10) 

p g = * + (l)o(l) ... <r(5)a(5)P o (6)o(6)P o (7)p(7)p + (8)a(8)6 + (9)p(9)P_(l0)Q(l0)P_(ll)p(ll) 

(r'(l2) p( 12) ... 6 + '(l6) p( 16 ) (6. 11) 

p 9 = 6 + ( 1) a( 1) . . <r(5) a(5) P q (6) q(6) P q ( 7) p(7) 6 + '(8) a(8) P + (9) p(9) P_( 10)«(10) P_( 11) p( 11) 

er'(l2)p(l2). . . 6 + '(l6)p(l6) (6.12) 

p 10 = 6 + (l)a(l) ... »(5)a(5)P o (6)Q(6)<r(7)p(7)P + (8)a(8)P + (9)p(9)P_(l0)a(10)P_(ll)p(ll) 
<r'(l2)p(l2) ... 6 + '(l6)p(l6) (6.13) 

P U » 6 + (l)a(l) ... «K5)a(5) o*(6) a(6) P q ( 7) p(7) P + (8) a(8) P + (9) p(9) P_( 10) o(l0) P_( 11) p( 11) 
«r*( 12) p(l2) . . . 6 + '(l6)p(l6) (6.14) 
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The matrix elements of the energy for the S = 0 case are given as: 


H n = yUOpjIHIpj) (6.15) 

H l2 = Y (AOp 1 IH|p 2 + p 3 + p g + p 9 ) (6. 16) 

H 13 = \(AOp 1 lH|p 4 + p 5 + p 6 + p ? ) (6. 17) 

H u = y(AOpjlH|p 10 + p n ) (6. 18) 

H 22 = ^ AO { p 2 + p 3 + p 8 + p 9} |H| { p 2 + p 3 + p 8 + p 9} * ( 6 - l9 ) 

H 23 = ^ AO { p 2 + p 3 + p 8 + P 9} IHI { p 4 + P 5 + P& + P7} > (6. 20) 

H 24 = * AO {p 2 + P^ + P8 + p 9} lHl { p 10 + p l 1} > (*-21) 

H 33 = y(AO [p 4 + P 5 + p 6 + p 7 } IHI |p 4 + p 5 + p 6 + p 7 } ) (6. 22) 

H 34 = y(AO |p 4 + P 5 + P 6 + P 7 } IHI [p 10 + p u }) (6.23) 

H 44 = y(AO {p 10 + p n } IHI {p 10 + p u }) (6.24) 

Y = y/Z 


As an illustration of how the spin operator method works the element 
y(AOp 2 /H/p^) will be worked out in detail. All of the one-electron orbitals are ortho- 
normal and it has been assumed that there is no overlap between the orbitals localized 
about the two Mn ions. The spin operator used can be expressed as 


. N w 

0 = ^ {l + c ! I (S' ) + c I 

Y 1 1 . k_ J A :.>l 


j A »k A «b>^ 


S L>♦•••} 



Since O is to operate on a spin product function, we may speak of the electrons as 
being distinguishable. Electrons 1, 3, 4, 5, and 10 comprise set A. Electrons 12, 
13, 14, 15, and 16 make up set B. N is the number of electrons in set A or B and 
M is the number of p spins assigned to the electrons in set A or the number of a 
spins assigned to set B. Therefore, the term 
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N 

I s: s* 

V B >A B 


in O means that runs over the electrons in set A and kg is to run over the elec¬ 
trons in set B. (AOp 2 /H/pg) written out is 


AO {& + ( 1) a( 1) 6 (2) a(2) w + (3) a(3) w_(4) a(4) <r(5) a(5) P q ( 6) a(6) P q ( 7) p(7) P + (8) a(8) 

P + (9) P(9) P_( 10) a( 10) 6_( 11) p( 11) ff'( 12) p( 12) »_'(13) p( 13) x + '( 14) p( 14) 6_'( 15) 
P(l5)6 + '(l6)p(l6)} I HI {6 + (l) a (l)6_(2)Q(2)w + (3)a(3)*_(4)a(4)a(5)a(5)P o (6) (6.26) 

a(6) P q ( 7) p(7) P + (8) a(8) P + (9) p(9) w_'( 10)a( 10) P_( 11) p( 11) <r'( 12) p( 12) *.'( 13) 
P(l3)w + '(14)P(14)6_'(15) P(15)6 + '(16)P(16)} 


where the orbitals which do not match are underlined. 

The 1 from 0 combined with the 1 from A gives for a contribution 

(P_( 10) 6.(11) | w_'(10) P_(ll)) 

Due to the spatial orthogonality of the orbitals and due to spin orthogonality, the only 
contributing permutation from A that can be combined with the 1 from 0 is the 
identity permutation. 

The term 



combined with the 1 from A makes no contribution due to spin orthogonality. The 
only term to give a non-zero result is Cj Sj Q s| 3 which must be combined with P^ 0 13 
P 10 i 3 from A > the combination yielding 

-0^6 (11) P_(13) IP (11) w_'(l3)) 

There will be no other terms entering into this matrix element. Therefore, 

(AOp 2 )IHIp 5 ) = (1 - Cj)(6.( 1) P_(2)|P_( 1) w.'(2) (6.27) 
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and Cj is - 1/5. By proceeding in exactly the same manner, all of the matrix elements 
of H for the S = 0 case can be set up. 

Let the Hamiltonian matrix for the singlet states be denoted by H Q and the 


H = H. + d‘ 4 > 


The matrix D has the following form 


O O O O 


O D,, D,, O 


Now suppose >|» 0 is the eigenvector of H Q corresponding to the lowest eigen¬ 
value of H Q which we shall call E Q . Then 

(4. 0 IHl+ 0 ) = (4* 0 IH 5 14» 0 ) + (+ 0 lDl+ 0 ) (6. 


If is the eigenvector of corresponding to the lowest eigenvalue E^, then by the 
variational principle 


(* 5 IH 5 I* 5 ) $ (^qIh 5 I4» 0 ) 


E 0 > E, + D 
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lyfwrli *'t. 'I. 4 . n. r « w . * respectively, and ; tor k = 12 . 

II. "ifT)§r m **>. -&r*. *1 a '. <**>4 «,' respectively. O33 is oblaiaed t>y replacing 

~& t i^ 1 i cv«'r«v?Sft-» , c 4 n 0 ^ except when 6 + ami 5 ^' occur as the restilt 

+ (6ji.)P 0 (2)iP 0 n)^„'(2)) * u_(i)p + (2)ip + ( 0 v(iy 

+■' g- ( 6 _ 1 (1) P _ (2) | P _ {1) » J 2)) + U '< 1) P q ( 2)|P J l| tr_(2)) + <6.'(l) P f (2)|P f ( 1) ir_(2)) 

- + |(6 r d)P t (2)|P.{l)ir + *(2)) + (6 . (1) P 0 (2)|P o < 1) *.'U)") i (6.(1) P_U)| P_( 1) tr + '(2)) 

+ |(6 + , (l)P + (2;)P + U)Tr + (2)) + (6 + '( l) P 0 (2)|P„(1) tr + |2|) f (6 + J(l) P,-(2)|P-_(l) ir + (2)) 

: (6. 34) 

Ab .5. 

°44 = ? £ (4>/(i)P 0 (k)|P (1)+ k J (k) + 't 1 (*. (1)P (k)|p u)4> k (k)) 

- 5 k = 12 - K J ° K 5 k = l - ° ° K 

+ -yCHl)P fi (2)|P o 1l)aH2)) + 2(tr(l)P + (2)|P + (lJ 0 '(2)) * 2(u( l) P_(2) |P_( l) <r l U)) 

( 6 . 35 ) 

Tfi'trprincipal parts of and are exchange integrals between the oxygen 
■ orbitals and the Mn orbitals; these integrals are all. positive. The other integrals en¬ 
tering D.^ , and are three-center integrals which cannot be of major import as they 
involve the orbitals of both Mn ions. Although these three-center integrals -are by. no 
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a 16 a 5 

D 22 = I 2 (« k '(l)P.(k). , P.(l)« k '(k)) l |^Z^(4' k (l)P.(k)|P_(l)<j> k (k)) 

16 5 

♦ l k ^ 12 ( *k’ (l)P + (k),P + (l) «k ,(k)) + | i (<*» k (D P + (k)IP + (l)4» k (k)) 

+ H(p_(i)6_U)|6_'(l)P_U)) + 2(6_(l)P o U)lP o (l)6_'(2)) + 2(6_(l) P + (2)|P + (l) 6_'U)) 

+ ^(P + (l)6 + (2)|6 + ’(l)P + (2)) + 2(6 + (1)P 0 (2 )|P o (1)6 + '(2)) + 2(6 + ( 1) P_(2)|P_( 1) 6 + '(2)) 

(6.33) 


We mean by $ k for k = 1.2, 3, 4, 5; 6 + , 6 . w + , w , a respectively, and for k = 12, 
13, 14, 15, 16; <r', w_\ w + \ & ', and & + ' respectively. D 33 is obtained by replacing 
6 + by » + and 6 + ' by w + ' everywhere in D 22 except when 6 + and 6 + ' occur as the result 
of a summation. 

°23 = |(Ml)P_U)lP_(l)w_’(2)) + (6_(l)P o (2)|P o (l)v_'(2)) + (6_(l) P + (2)|P + ( 1) *_'(2)) 

+ |(6_'(l)P.(2)|P_(l)w_(2)) + (6.’(1)P 0 (2 )|P o ( 1) W _(2)) + (&_’(l) P + (2)|P + (l) tt_(2)) 

+ |(6 + (l)P + (2)|P + (l)w + ’(2)) +(6 + (l)P 0 (2)|P o (l)w + '(2)) + (6 + (l) P_U)|P_(l) w + 'U)) 

+ |(6 + *(1)P + (2)|P + (1)w + ( 2)) +(6 + ’(l)P 0 (2)|P 0 (l)n + (2)) + (6 + '( 1) P_(2)|P_(l) w + (2)) 

(6. 34) 

D « = 5 kl 5: i2 ( *k' (,,P o (k > IP o< 1 >*k , "‘l * | k | i ( V ,,p . , k » lp 0 ( , )* k ( k » 

+ -^-(<Kl)P 0 (2)|P 0 (l)^'(2)) + 2(«r(l)p + (2)|p + (l)«r'(2)) + 2(a( 1) P_(2)|P_( 1) «r'(2)) 

(6. 35) 

The principal parts of D 22 and D 33 are exchange integrals between the oxygen 
orbitals and the Mn orbitals; these integrals are all positive. The other integrals en¬ 
tering D 22 and D 33 are three-center integrals which cannot be of major import as they 
involve the orbitals of both Mn ions. Although these three-center integrals are by no 
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means zero, they rre overshadowed by the exchange integrals. Therefore D 22 and 
D 33 as well as D 44 are positive. Since D 44 is positive, if D is to have a negative 
eigenvalue it must come from the 2 x 2 secular equation involving D, 2 , D 33 , and D 23> 
The secular polynomial coming from this part of D is 

^ * (D 22 + °33> K 4 D 22 D 33 ' ( D 23^ = 0 ( 6 - 36) 

We can only have a negative root if D 2 2 D 33 < ^°23^ since the sum of the roots is 
positive. By means of the same arguments which established the positive character 
of D 22 and D J3 we can see that D 23 is less than either D 22 or D 33 - Therefore, the D 
matrix has no negative roots and consequently the energy of the lowest lying singlet 
state as found by this configuration interaction cannot lie below the lowest S = 5 state. 

Up to this point configurations in which an electron on the oxygen ion goes into 
a 4s orbital on one of the Mn ions have not been included. The 4s orbitals are more 
extended than the 3d's and now we would have to admit exchange integrals between a 
4s orbital localized about one Mn ion and the 3d orbitals localized about the other Mn. 

By examining the character of D 22 and D 33 it appears that the inclusion of these con¬ 
figurations could not lead to a lower lying singlet energy than that of the S = 5 states. 

There is one quite definite conclusion which can be drawn from this investiga¬ 
tion. Within the framework of the approximations made here It seems clear that when 
an electron is removed from the O' ion. leaving an unpaired spin on the oxygen, that 
this unpaired spin couples the non-overlapping charge densities on the Mn ions ferro- 
magnetically. In view of this result, it appears likely that some other appreaen would 
be more profitable as an explanation of antiferromagnetism. Of course, it may well 
be that the singlet state considered here in no way resembles the antiferromagnetic 
state. 

Before abandoning this simplified model there is one feature of the configura¬ 
tion interaction treatment that bears examination. That is the fact that this approach 
has placed all of the emphasis on the probability amplitudes of the various configura¬ 
tions, the energy being stationary with respect to the variation of these parameters. 

It is quite possible that the one-electron wave functions are really quite different for 
the singlet and S = 5 states. Furthermore, the use of orthogonal orbitals may force 
one to carry out a rather extensive configuration interaction to accurately describe the 
correct state of affairs. An investigation is being carried out at present using a mix¬ 
ture of orthogonal and non-orthogonal orbitals but considering only the ground configura¬ 
tion. The energy will be stationary with respect to variations of parameters in the one- 
electron functions and these functions will be determined separately for the states of 
different multiplicity. Preliminary results obtained for this approach have been en¬ 
couraging. 
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7. AN EXTENSION OF HUND 'S RULE 


In the course of configuration interaction calculations with orthogonal orbitals 
it is sometimes interesting to know what the multiplicity is of the lowest state of some 
given class of states. In general it is impossible to answer this question without cal¬ 
culations. Under one set of circumstances it is possible to answer this question with¬ 
out any calculation. We know of Hund's rule of atomic spectra. Part of this rule states 
that of all the levels arising from a given spacial configuration the state of highest mul¬ 
tiplicity lies lowest. A rule similar to this can be proved in general. 

We imagine n distinct orthogonal orbitals Uj, u^, • • • u n - We put one elec¬ 
tron in each of these orbitals with either spin up or spin down. This gives rise to 2 n 
states. These states can be formed into states of definite multiplicity ranging from 
singlets (for n even) to states of multiplicity n + 1. (The presence of other orbitals 
which are doubly filled does not influence the validity of the following arguments and 
can be dropped from considei ition. ) We shall now show that there is no state that lies 
lower than the state of highest multiplicity. We lose no generality by restricting our¬ 
selves to all those states with the same z component of total spin. Of these only one 
corresponds to the state of highest multiplicity whereas for any other multiplicity 
there is in general more than one state of this multiplicity. Since states of different 
multiplicities cannot interact we are confronted with the problem of diagonalizing for 
each multiplicity the matrix of the Hamiltonian 

H = £ - V 2 ♦ f(i) ♦ I g u (7. 1) 

i>j 

Here f(i) is the one-electron potential and g. . is the interelectronic Coulomb interac- 

(i \ 

tion. From the Dirac vector model we know that for a given multiplicity the matrix 
of the Hamiltonian can be written down very simply. 

A, = co M ,. + I - J (P ) w (7.2) 

ij 

Here represents the Hamiltonian matrix and (Py)^ are unitary matrices which 

represent the permutation of the spatial coordinates of electron i with those of elec¬ 
tron j. Jy is the customary exchange interaction = Uj *(l) g u u t (2) u (2) 

dTj dT^ which is greater than zero since it represents the interaction of a charge den¬ 
sity with itself. We first consider the state of highest multiplicity. This has a Ham¬ 
iltonian matrix which has only one row and one column. The matrices representing 
the permutations are matrices which have one row and one column. The matrix 
element in all these matrices is unity. Therefore the energy of the state of highest 
multiplicity is (apart from the constant in (7. 2) which is the same for all states) 
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E n ♦ 1 = * Z J ij 

For any other multiplicity the matrices are of larger size. Let us consider the 

diagonal energies of all the states of a given multiplicity. These energies are the di¬ 
agonal elements of (7. 2) 


l- W., 


E „*. *- I'n I- 


We have therefore proved that there is no state lower in energy than the state of highest 
multiplicity. One might ask if going through the process of diagonalizing the Hamilton¬ 
ian matrix (7. 2) will push some state lower than the state of highest multiplicity. This 
is not the case since in diagonalizing the matrix (7. 2) we use a unitary transformation 
which does not destroy the unitary property of (P^)^. This unitary property was all 
that was necessary to prove the theorem. In general, in order to get some state other 
than the state of highest multiplicity lying lowest it is necessary to take states into ac¬ 
count which occupy the orbitals in a different manner possibly filling some of our n or¬ 
bitals twice and leaving others empty. 
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3. ON THE O z GROUND ELECTRONIC WAVE FUNCTION 

Microwave absorption measurements on O 1 ^ O*^ gas^ exhibit a hyperfine 
structure. The observed spectrum agrees with that predicted on the basis of a perturb¬ 
ing Hamiltonian of the form bT • £f ♦ cl S ; <f is the electron spin angular momentum 
17 z z 

of the molecule, and 1 is the O nuclear spin (I = 5/2). According to a i'neory by 
Frosch and Foley, ^ 2 ' experimental values of 


<P z /r l > = /( 3 cosZ 3 x ~ 1 ) p(r) dr 


( 8 . 1 ) 


and p(0) can be deduced from those of the parameters b and c. The angle between the 
internuclear axis (the z-axis) and the radius vector ? from the O 1 ^ nucleus to an elec¬ 
tron is denoted by x. and P 2 (x) = (l/2)(3x 2 - 1). <P 2 /r 3 > is the value of P 2 /r 3 aver¬ 
aged with respect to the number density p(r) of electrons of unpaired spin of the mole¬ 
cule, and p(0) is the value of p at the O 17 nucleus. In the following we compare the 
experimental values - 1. 32 and 0. 187 (in atomic units) of <P 2 /r 3 > and p(0), respectively. 

obtained by Miller and Townes^' with values determined by Meckler's O, normal 
(31 Z 

electronic wave function. ' 

<P 2 /r 3 >may be defined more precisely as 


<Pj/r 3 > = /*•(? 


P 2 (cos Xj), 


$ dr 


( 8 . 2 ) 


1 electron measure in units of )(. The 
summation is over the N electrons of the molecule, and the integration is ever the co¬ 
ordinates of all electrons (and is to include summation over spin coordinates). $ is a 
normalized electronic wave function, the solution of Schr5dinger's equation for 0 2 with 
a spin independent Hamiltonian and with fixed nuclei. $ is an eigenfunction of if 2 and 
S z with eigenvalues S(S + 1) and Mg = S; jf = XjS- Upon integration over the coordinates 
of all electrons but one, and summation over all spin coordinates, (8. 2) reduces to (8. 1) 
in which p(r) is defined by 


p(r) = N JV(X, X 2 .X N ) 2s z $(X. X £ .X N ) d«r dX 2 . . . dX N (8. 3) 

and satisfies J p(r ) dr ~ f$+ 2S Z $ dr = 2S. In case $ is expressed as a linear com¬ 
bination of determinants of orthonormal orbitals, p(r) is a quadratic form in the orbitals, 
the cross-terms (overlap terms) arising from products of determinants differing by one 
orbital. 
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Quite generally p(r ) can be employed to find the expectation value of F: 

J}* F $ dT = J f(r) p(r) di\ where F is an operator of the type F = f(r^) 2 s j Z - 

E. g. . in the theory of neutron-electron magnetic scattering^ such an operator occurs 
with f(r) = 1 e'* 1 ' ' r where K is the neutron scattering vector. 

The normal O^ electronic level is a 3 ]T‘. Meckler^ approximates if for 

this level by a linear combination of nine symmetry states of 3 T symmetry, each a 

“g 

linear combination of determinants of orthonormal orbitals, these orbitals being sym¬ 
metry orbitals (LCAO molecular orbitals) from one configuration. We write Meckier’s 
if as 


*M “ £ C u 


(8.4) 


where, in terms of the determinants specified in Table 8-1, the symmetry states are 
given by 


<t» a = j(A + B - C - D) 

4> b = i(B + C - A - D) 

$ g = 2* 3 / 2 (E + G - F - H + J + K- I-L) (8.5) 

♦ h = 2* 3 / 2 (F + G- E- H + I + K- J-L) 

4>. = 2" 3 / 2 (E + F - G - H + I + J - K - L) 

and the one-determinant symmetry states A , $ ,, 4> , and The coefficients c as a 

/jxCcie i 

function of R are given by Meckler. The result of inserting (8. 4) into (8. 3) is, in 
terms of Meckler : s symmetry orbitals (see Table 8-1), 

p(r) = C,|* o | 2 + C 2 |x o | 2 ♦ C 3 |x t | 2 ♦ C 4 I4* ± | 2 (8.6) 

where 

Cj = a + p+ Y- 6 
C 2 « a + p * Y + 6 

C 3 = lc a+ c b | 2 + 2lc c l 2 + 2|c e ! 2 + a - p + y + 6 
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Table 8-1 

Determinants of the symmetry states* 


A 

B 

C 

D 

E 

F 

G 

H 

I 

J 

K 

L 

♦c 

*d 

*e 

♦f 


2 p«r 


2 p? 

x o 



*P" U + 2pir u 

4 > + 


*The normalized determinants A, . . . , with Me = 1 are defined by the symmetry or¬ 
bitals on the principal diagonal; these consist of filled Is? and 2s? orbitals accounting 
for eight of the sixteen electrons, the remaining orbitals being specified by the above 
table, a + or - entry indicating that the space orbital is occupied by an electron with 
+ or - spin. The space orbitals are labeled both in standard notation and in Meckler's 


notation. Meckler lists the nine 3 £’ 
eigenvalue Mg = 0; these yield the 


symmetry states as eigenfunctions of S z with 
when operated on by (S x + iS^)/ ->/T. 
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\ + 6 


° = 4 K g ' C h + C i' P = 4 1 * c g * C h + c i' 


4^g h "i 


c ,l 6 = T l c n ♦ c h + c il 


No overlap terms appear in p since all determinants involved differ by at least two or¬ 
bitals. This follows from the fact that no two space orbitals available for occupation 
outside the lscr and 2s<r filled shells are of the same symmetry. Values of the C's and 
their component parts are given in Table 8-2. At normal internuclear distance R = 2. 28 


Table 8-2 

Values of the C's of Eq. (8. 6) and their component parts for Meckler'r 
0 2 normal 3 £’ electronic wave function. Note that Z k C k = 2.00. 


R 

2 . 0 

2. 5 

4. 0 

00 

C 1 

- 0.0018 

- 0.0036 

0.4424 

1/2 

C 2 

0.0024 

0 . 0066 

0. 4425 

1/2 

C 3 

1. 9632 

1.9183 

0. 5960 

1/2 

C 4 

0. 0361 

0.0788 

0.5181 

1/2 

a 

0.0003 

0.0013 

0.2226 

1/4 

P 

0.00002 

0.0002 

0.2202 

1/4 

Y 

0.0060 

0.0113 

0.00004 

0 

6 

0 . 0081 

0.0164 

0.00008 

0 

21 c | 2 
c 

1.9400 

1.8694 

0. 3804 

1/4 

21 c/ 

0.0216 

0. 0479 

0.3208 

1/4 


0.0084 

0.0183 

0.2129 

1/4 

2'%| 2 

0.0007 

0. 0041 

0.1998 

1/4 

|C a - c/ 

0.00000 

0. 00004 

0.00001 

0 

|c a * c l/ 

0.0005 

0 . 0018 

0 . 0001 

0 
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the symmetry state * . which corresponds to the molecular orbital solution for the 
c 2 

normal 0 2 molecule, predominates (l - Ic^l « 1). The p for which involves only 
2 pWg orbitals x+. is not much different from the p for 

The value of <P 2 /r 3 > calculated from p of at equilibrium internuclear 
separation is only about 27 percent of the experimental value - 1. 32. When Meckler's 
orbitals in p are replaced by Hartree-Fock orbitals, the value of<P 2 /r 3 > is about 90 
percent of the experimental value. These values of <P 2 /r 3 > are discussed in more 
detail in the next paragraph. It seems unlikely that if Hartree-Fock or other reasonable 
types of atomic orbitals were used in Meckler's calculation in place of his Gaussian- 
type orbitals, the C's of (8.6) would be significantly different. The value of <P 2 /r 3 > 
will then depend largely on the types of orbitals. One expects that calculated values of 
the molecular energy and <P 2 /r 3 > as a function of the atomic orbitals employed are 

not very well correlated, since <P,/r 3 > , but not the energy, is sensitive to values of 
L 17 

the orbitals primarily in the vicinity of the O nucleus. 

Let us consider first <P 2 /r 3 > calculated for We then have, in terms 

of the atomic orbitals p + about the O 17 nucleus and p + ' about the O 1 ^ nucleus 

<P 2 /r 3 > = 2 fP IX+! 2 dr = (2 • *£) A [l - c 2 (R) ♦ c 3 (R)] (8. 7) 


where 

A = /(P 2 /r 3 ) Ip +! 2 d7 

c 2 (R) = 2 f( P 2 /r 3 ) p t * p/d?/A (8. 7') 

« 3 (R) = /(P 2 /r 3 )|p t '| 2 dP/A 

and K(R) is a normalization constant. A, which is - <P,/r 3 > for an isolated O atom 
4 3 c 

in its ground level 2p P, is independent of R; it is negative, since the angular part of 
the integral is negative. For Meckler's Gaussian type p orbitals, A is - . 304 and, if 
we ignore the other factors in (8.7), yields a theoretical value of <P 2 /r^> only 23 per¬ 
cent of the experimental value. The theoretical value is increased somewhat by the 
remaining factors in (8. 7). The contributions i 2 and i 3 involving the O l6 orbital are 
small, and should decrease as R increases; in fact, for Meckler's orbitals and R = 2, 
c 2 = - 0. 03 and ICj/c^J ~ 0. 04, so that the factor (1 - t 2 + t 3 ) increases* the theoretical 

*Actually c 2 is the algebraic sum of three terms whose magnitudes for R = 2 are of the 
order ten times U 2 I, and it is not obvious that Ic^l decreases monotonically with in¬ 
creasing R. 
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value by about 3 percent. K 2 for Meci.Iei *s orbitals has the values 1. 25, 1. 09. and 
1. 00 at R = 2. 0, 2. 5 and 4. 0, respectively, so that the normalization of x + increases 
the theoretical <P£/r 3 > by about 15 percent at normal internuclear separation. The 
net result is thj:t <P 2 /r 3 > for the molecular orbital ground state of O^ computed us¬ 
ing Meckler's Gaussian type orbitals is still or.ly 27 percent of the experimental value. 

A for a Hartree-Fock O orbital^ la ^ is 71 percent of the experimental <P£/r 3 >. If, 
however, one includes the other factors of (8. 7), assuming these are unchanged in value, 
the discrepancy cam be reduced to about 10 - 15 percent, perhaps even less when these 
factors are properly computed. The value of <P£/r^> for is in slightly poorer 
agreement with experiment than that for $ . Now p involves integrals like (8. 7) with 
the 2p» u , 2p<j-g, and orbitals. Since the latter integrals contribute to 

with small coefficients, they need not be investigated with great precision. For 2pw 
the integral differs from (8. 7) in that K 2 is replaced by Meckler's L 2 < 1 and -* - c^, 
and so is smaller in magnitude than (8. 7). For 2po- orbitals the sign of the integral is 
reversed. A being replaced by - 2A. It is clear that <P 2 /r 3 >^ does not differ appre¬ 
ciably from <P 2 /r^>^ * 

Let us next compare the experimental and calculated values of p(0). For 
p(0) is ve.'y much less than the experimental value. This is not surprising, since only 
symmetry orbitals containing s orbitals have an appreciable density at the O 17 nucleus, 
and involves s orbitals outside the filled ls<r and 2sa shells only in the 2p«r orbitals, 
and then only in virtue of the orthogonalization of the 2ptr to the ls«r and Zs<r orbitals. 
Table 8-3 shows values of p(0) for some atomic orbitals. Values of + ls (0)| 2 and 

Table 8-3 

Values of p(0) for some atomic orbitals 

0 (5) ° 2 (3) 

Hydrogen-Like Hartree-Fock v ' Meckler^' 

Is 0. 318 Z 2 148.2 45.4 

2s 0.318 (Z/2) 2 7.604 0.363 

„ 2s (0)l 2 f° r Meckler's orbitals are smaller by factors 0. 306 and 0. 048, respectively, 
than the corresponding values for the Hartree-Fock O orbitals of Hartree, Hartree 
and Swirles. ^ Even when the latter replace Meckler's atomic orbitals in his symmetry 
orbitals, p(0) is less than 0.001 for smaller by a factor of about 100 than the ex¬ 
perimental value 0. 187. (We have neglected the reorthogonalization and renormaliza¬ 
tion that should be performed for these substituted atomic orbitals. ) 

Apparently the normal electronic state contains small admixtures of other 
symmetry states than those considered by Meckler, and these additional symmetry 
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states have large effects on p(0). The situation can be clarified qualitatively to some 
degree. A rough criterion for the magnitude of the amplitude c^ of a symmetry state 
is the sum of the energies of the component orbitals, the higher the energy, the smaller 
Ic |. Only the <r orbitals of a symmetry state contribute appreciably to the density of 

17 

unpaired spin in the vicinity of the O nucleus, since only these contain s orbitals. 
Thus, only symmetry states with unpaired a orbitals can contribute effectively to p(0). 
We can get some idea of the relative effects on p(0) of various ns* orbitals by noting 
that for hydrogen-like atomic orbitals 41^, * >4< n Q(0)/^*i Q (0)l 2 = n" 

On the basis of these ideas we may expect that only symmetry states with low 
energy containing unpaired nstr orbitals with small n give important contributions to 
p(0). Examples of such are symmetry states which differ from one of the in that a 
1 s<r or 2s a orbital is replaced by an n la of low n * 2 and I / 0, or in that a 2p* is re¬ 
placed by an nstr with small n > 3. (Note that if the excited symmetry state is to have 
the same symmetry as a ^ and involve determinants which differ by only one orbital 
from thos of the the excited orbital must have the same symmetry as the unexcited. 
Also, a cross term in p(r) vanishes at r = 0, unless the orbitals which overlap are 
both s a ) 

Approximate *ave functions in molecular theory are commonly determined 
by a variational procedure which makes the energy stationary. The result is that al¬ 
though the energy is insensitive to small variations in the wave function, other quanti¬ 
ties, such as <P 2 /r*> and p(0), may be relatively quite sensitive. One should then 
not expect amounts of mixing of symmetry states which yield energies in satisfactory 
agreement with experiment to necessarily yield other molecular quantities in such satis¬ 
factory agreement; p(0) is an example of such a quantity. On the same basis, types of 
atomic orbitals which lead to reasonable energies may lead to poor values of other mo¬ 
lecular quantities; both <P 2 /r^> and p(0) are examples of such quantities. 
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9. THE ETHANE MOLECULE 


As indicated in the last Report, the severe restriction on the amount of con¬ 
figuration interaction to be considered necessitates an optimization of the molecular 
orbitals, a best choice of linear combinations of atomic orbitals. This problem is faced 
also by Dr. Kaplan in his treatment of the NH^ molecule where even though quite a bit 
of configuration interaction is handled, the amount is still not enough to permit complete 
equivalence of all linear combinations of atomic orbitals. Roothaan^ has given a pro¬ 
cedure for determining best LCAO's for the case of a closed shell state, a single de¬ 
terminant many-electron wave function. His equations are the algebraic equivalent of 
the Hartree-Fock differential equations and are similarly characterized by the dis¬ 
couraging prospect of a self-consistent calculation. However, there is Whirlwind -- 
the high speed computer now of large capacity -- and Roothaan's procedure can be mech¬ 
anized and performed by the machine so that it becomes a fast, repeatable routine. We 
wish to present flow chart of the machine program now being developed but first we 
will re-express Roothaan's equations in a slightly modified form: 

To begin with, there are the one-electron functions which are to be combined 
in the best possible ways. Let them be orthonormalized. Here we differ with Roothaan 
who did not require orthogonality. The lack of orthogonality leads to a more compli¬ 
cated type of determinantal equation to be solved and we see no point in adding this diffi¬ 
culty to the non-escapable one of self-consistency. We begin with a well-considered 
first guess to the molecula/ orbitals, an orthonormalized set. A general linear com¬ 
bination of these functions is still a general linear combination of the original atomic 
functions. 

This original set of functions is denoted as v^ and the number in the set by N. 

If there are 2n electrons in the system, we are to find n linear combinations of the v^'s 

♦i * I % V 

such that an antisymmetrized many-electron wave function formed with each doubly 
occupied will be associated with a minimum value of the energy. The expression for 
the energy is 


where 


E * 2 I v H v„ + I v G Ow 

M-v pvXff 

V ■ /V (1) [• v * v(,) ] v v (1) d3x > 
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G W = 2 ^ wlXa ) " 

(nwlXff) = f v *(l) v (1) -i- v *(2) v (2) d 3 x. d 3 x ? 

J \k v r i 2 X o- 12 

r* * 

P = > c. c. 

r (lV t-‘ 1(1 IV 

i 

or in matrix notation 

p = C + C (+ = Hermitian adjoint) 

and the orthonormalization of the ^'s is implied by 
CC + = 1 

C is an n x N rectangular matrix, p is a square N x N. The n rows of C, 
which are the complex conjugates of the n columns of C, are the n lowest eigenvectors 
of the matrix 


That is 


K. = H. ♦ V p G^ v 
X cr Xo- L r pv X cr 

pv 


KC + = AC + A diagonal 


We are to pick a p, form K, diagonalize it, selects its n lowest eigenvectors, form 
p, re-cycle. The thing to realize is that H and G do not change from cycle to cycle. 
They are the input to the machine, the data to be used over and over again. The flow 
chart is illustrated on p. 48. 

The machine is to do matrix-vector multiplication (symbolized as pG), matrix 
addition, matrix diagonalization, magnitude selection and comparison, and many cycles. 
Even with the type of factorization by symmetry discussed by Roothaan there is a lot to 
be stored in the machine, but the addition of more electrostatic storage tubes and a 
magnetic drum to Whirlwind has made for the pleasant atmosphere of a buyer's market 
in storage space. 

The matrix diagonalization procedure, which is just a sub-routine in the whole 
self-consistent process, has been programmed and will exist as a separate routine in 
the Whirlwind library. The present program can handle up to a 31 x 31 real symmetric 
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matrix and a 15 x 15 complex Hermitian matrix. The rest of the L.C.A.O. determina¬ 
tion is now being programmed and will be first applied to the NHj molecule. 

Alvin Meckler 
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10. LIMITED CONFIGURATION INTERACTION TREATMENT OF THE NH, MOLECULE 


The calculation of matrix elements indicated in the last Progress Report is 
still underway. The matrix elements have been expressed in terms of "directed" mo¬ 
lecular orbitals. Due to the symmetry properties of these one-electron functions, 
there are only thirty-five independent two-electron integrals and five independent one- 
electron integrals. Each of these integrals, however, is a combination of many similar 
type integrals involving atomic wave functions. At present the component atomic in¬ 
tegrals are being expressed in a form necessary for the evaluation of their radial sub¬ 
integrals by I. B. M. machine methods. 


H. Kaplan 
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11. THE WATER MOLECULE 


The tabulation of the a n 's mentioned in the last Quarterly Progress Report 
is nearly complete. We now have the following ra n 's; 

r = 0(. 1) 1.2(. 2) 4. 0(. 5) 7(l) 12 

a * 1.4, 1.5, 1.6, 1.8, 2.0, 2.2, 2.6, 3.0, 4. 0 and 5.0 

k = 1, 2 

n = 0, 1, 2, 3, 4, 5 (for k = 1) 

n = 0, 1, 2, 3 (for k = 2). 


The interelectronic interaction integrals involves terms of the form 


oo r i 

/ pab^i*[ViTi / «»cd (r 2 )r 2 ,dr 2 + r i I ^rnr dr i (ii i) 

o L r ! o rj r 2 J 

These can be written as 

/ l> ab< r l)[^Tn / P cd (r 2> dr l * / >'cd< r l> j^TTi / P ab < r 2> r 2 ,dr 2j dr , 


( 11 . 2 ) 


We shall call the quantity in square brackets in (11.2) the inner potential. 

We have written and tested a program for computing the inner potential, 
where p ab is the product of two radial probability amplitudes, for the I. B. M. Card 
Programmed Calculation. Our mesh is 


r * 0(. 0125) . 05(. 025) . 2(. 05) . 6(. 1) 1. 2(. 2) 4(. 4) 10 . 


The program takes three minutes of machine time per integration. The results of this 
indefinite integral are printed and punched on cards ready for use in definite integrals. 
This program is now in use. 


G. F. Koster 


12. CONFIGURATION INTERACTION APPLIED TO THE HYDROGEN MOLECULE 


Calculation of the several matrix elements discussed in the previous Progress 
Report is continuing. The odd one-electron function obtained from the Hartree equa¬ 
tion containing an "average field" due to an electron in the even orbital has been deter¬ 
mined. The kinetic energy and nuclear framework integrals appearing in the main di¬ 
agonal element in which the electrons'are in the-odd orbita-l-have--bees.iompu.ted and . 
there remains only the electron-electron interaction integrals to be found. 

E. Callen 
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13. CONFIGURATION INTERACTION FOR THE FLUORINE MOLECULE* 1 ^ 

In order to shorten the computation in the problem mentioned in the last 
Progress Report, it has been decided to use Slater-type atomic orbitals^ instead of 
Hartree-Fock atomic orbitals. Work on the problem is proceeding along these lines. 
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14. NUCLEAR ELECTRIC QUADRUPOLE INTERACTION IN THE KC1 MOLECULE 


The introduction of a polarization effect by appropriate distortion of atomic 
orbitals has been studied along lines indicated in the previous Progress Report. ^ A 
special case of the proposed form: 

* - O 1 * r> ' <> < 14 -» 

has been investigated quantitatively. In particular, we have chosen n = 1, i = 0, m = 

0, I' = 1. m' = 0. and in addition we have provided for normalization and a "mixing 
coefficient", C. Thus our wave function is: 

<J< * Ne" r (l + Cr v cos 0) (14.2) 


where N = normalization constant and C, v = variation parameters. (Unless otherwise 
stated, atomic units, a. u. , are used throughout; a Q = unit of length = . 529 x 10 8 cm, 
e 2 /2a Q * unit of energy = 13. 53 e. v.) 

Using (14. 2) we have considered two types of perturbing fields: (l) a weak 
uniform field and (2) the field of a point charge a distance R = 1. 4 a. u. away from the 
center of the unperturbed function. 

Case (l) 



z. 


The Hamiltonian is: 




H « H Q (r) ♦ H'(r, * ) = - V 2 - | + £ r cos 0 


and H'(r, £ ) is considered to be small compared to H q . We form: 


E 


h* 

y 


H * dr 
4< dT 


(14. 3) 


Our problem is really that of the second order Stark effect for the normal hydrogen 
atom — which has been solved exactly^ yielding a polarizability: a » 2. 25 a. u. (4. 50 a Q ^ 
in any other system of units). Using (14. 2) and varying C and v we obtain a minimum 
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of (14. 3) for C = -.724, v = 1.50, giving a = 2. 25 a. u. --the correct value. For 
comparison we have used: 


* n i hs * c i*j P z] 



*2pz 


r cos e - r/2 
4 ^/zii 


(14. 4) 


Minimizing (14. 3) with respect to Cj we obtain a c 1. 48 a. u. (2. 96 a Q 3 ) — a value in 
error by 35 percent. The resulting expressions for and 4* are given below for Case 
(l). (We have separately normalized the angular dependent term of (14.2) so that the 
coefficient of the distortion terms are directly comparable!: 

+1 * N 1 [*L - • 9,5 '+2pz] N, - 1 - .««/* 

+ - N[e' r - 1.45*( r1 ’ 5 °J2IL? e N = .564 - . 369* 2 

It should be noted that for Case (l) (14. 2) is just of the form: 

4 > = +°(1 ♦ AH ' 4 . ..) 

which was suggested by Lennard-Jones^ on the basis of an approximate treatment of 
second order perturbation theory. Also Slater and Kirkwood^ 4 ' have treated the normal 
hydrogen atom in a weak uniform field using a function of the form: 

4 * « 4»°0 ♦ H'R(r)) 

They obtained a differential equation for R(r) su’d their resulting wave function was: 

4» c N [e' r - 1/ |(r + y-) cos 8 e' r J] N - .564 - . 282* 2 
and a « 2. 25 (4. 50 a Q 3 ) - the exact value. 
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Case (2) 


f?' = S I (i) k P k (cos ») r < R 

p, ■ p £ ( F )k P k (cos 0 ) r > R 
k 

Here we have a strong field case with H'(r, R) the same order of magnitude as H Q (r). 
Again using (14. 2) and minimizing the energy with respect to C and v for R = 1. 40 we 
find C » . 385, v = .65. For (14.4) we find Cj = . 109 and the respective functions are: 


The Hamiltonian is: 

H = H o (r) + H'(?\ R) « - v 2 - | | 



+ = N [e" r ♦ ■ 298 ( r c °^_^ e " r VI 

L yr6oo /J 

+1 * N l[+l S * 109 +2pz] 


N = .508 


. 989 


The exact solution of Schrodinger's equation for this potential is the hydrogen molecule 
ion with an energy at R = 1.40 of - 1. 15 a. u. + gives an energy - . 916 while <J/, gives 
- . 809, the first representing a 20 percent greater value than the exact value for H^ 
and the second a 30 percent greater value. However we are trying to solve a polariza¬ 
tion problem and not the H^ + problem and thus we have not brought in the symmetry of 
the positive charges. Our calculation only shows that the wave function of (14. 2) is 
more satisfactory than the case of (14.4) where we admix some 2p z function with the Is. 
The method of Slater and Kirdwood depends for the separability of their differential 
equation on linearity of H' in rectangular coordinates and spherical symmetry of + 0 
and appears not to be readily manageable for complicated cases. 

As a further test of the function (14. 2) the determinantal wave function for F~ 
is being set up using one-electron functions of the form (14. 2) and determining C and 
v by minimizing the total energy. 
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15. INTERSTITIAL DIFFUSION IN THE DIAMOND STRUCTURE 


We wish to derive the diffusion coefficient for interstitial diffusion ir the dia¬ 
mond structure by considering the elementary process of a transition of ! 1 t interstitial 
ion from its original site to a neighboring one. We here consider only inters. : tial diffu¬ 
sion, for the case of infinite dilution of the diffusing substance, constant temperature, 
and absence of electric fields. 

The diamond structure can be viewed as two interpenetrating face-centered 
cubic lattices, one having a parallel displacement of (a/4, a/4, a/4) with respect to the 
other, a being the length of the side of the unit cube The interstitial positions are 
found to have the diamond structure also, displaced by a body-centering translation 
from the original lattice. That is, if the carbon atoms arc found at (000), F. C., and 
(a/4, a/4, a/4), F. C., then the interstitial sites are located at (a/2, a/2, a/2), F. C., 
and (3a/4, 3a/4, 3a/4), F. C. The interstitial positions located with respect to it as 
as origin at n^ = (a/4, a/4, a/4), n 2 = (a/4, - a/4, - a/4), n^ = (- a/4, a/4, - a/4), 
n 4 = (- a/4, - a/4, a/4), and ihe interstitial position (3a/4, 3a/4, 3a/4) has nearest- 
neighbor sites located with respect to it by the negatives of these displacements. 

We assume that the transition probability is X. for a jump from one site to a 
nearest-neighbor site and O for any other jump. Then, if P(r y t) be the probability 
that site r^ is occupied by an interstitial ion at time t, we have 

4 

P(r t) = - 4XP(r t) + £ XP(r t n r t) 


where the upper sign holds if r^ is on the lattice containing (a/2, a/2, a/2) and the lower 
sign holds for the other case. A Taylor's series expansion provides immedi-tely the 
diffusion equation 

— (P, t) = ^ V 2 P(r, t) . 

8t 8 

The diffusion constant is thus D = X.a 2 /8. 

The neglected terms in the equation for either lattice of ions are the third- 
derivative terms, but after averaging over both lattices the only terms which survive 
are the fourth and higher derivatives. 

The assumption of infinite dilution of the diffusing substance can be relaxed 
somewhat, for it is found that if the probability of occupancy of target sites is considered 
the same diffusion equation results. The jumping probability X. will be a function of 
local concentration of diffusing substance, however, so the requirement of fairly high 
dilution would still hold. 


H. C. Schweinler 
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16 . A STLDY OF 2Z o IN ATOMS 


The results of calculations of the effective nuclear charge for potential in 
an atom have been submitted to the Physical Review for publication. 

A. J. Freeman 
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17. ATOMIC MULTIPLET SEPARATIONS 


Atomic multiplet separations have been computed for atoms with uncompleted 
3p and 3d shells. Simple analytic radial parts of the one-electron wave functions were 
assumed: 


V ... 2 -V 

- Z.r 


R, was obtained by taking a linear combination of single term expressions for the 2p 
■*P (i \ 

and 3p radial functions. The single term expressions are given by Professor Slater' ' 

as: 

r n* - 1 e - \z - S)/n*] r 


where Z is the atomic number and n* and S are given by simple rules. C g and are 
determined from orthogonality conditions with the 2p radial function and by normaliza¬ 
tion. If one assumes the radial parts of the one-electron functions do not vary from 
multiplet to multiplet for a given configuration, a multiplet separation is given by a 
sum of Coulomb and exchange integrals. This sum may be expressed in terms of the 
F k (n, 1; n\ V) and G k (n\ 1; n\ V) integrals. 

If orthonormal functions of the form, 




A -V 


are used, the integrals may be written: 


n + n,, n„ + n„ 


k (nl; n-i*) = 2 l C s C t C u C v r L S 1 U v (Z # ♦ Z,. Z u ♦ Z y ) 
stuv 




_ , , \"i n j 8 ni + "i , 00,_ 

' ( ’ e.'Vj k 
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Table 17-1 



Multiplet 

Calculated 

Observed 


Separation 

Separation 

Separation 

Silicon I 

3 P - *d 

10, 285. 5 cm" 1 

6, 198. 66 cm" 

Op) 2 

*d - *S 

15. 428. 2 

9, 095.43 

Phosphorus I 

4 S - 2 d 

17, 179.0 

11, 369. 1 

(3p) 3 

2 d - 2 P 

11, 452. 6 

7, 366. 1 

Potassium V 

4 S - 2 D 

20.914.6 

24, 118.5 

(3p) 3 

4 S - 2 d 

34, 857. 8 

39, 904. 5 

Titanium I 

3 F - *D 

6, 131. 02 

7, 069.62 

(3d) 2 (4s) 2 

3 f - 3 p 

7, 099. 02 

8, 324.8 


3 f - *g 

9,443.88 

11, 932. 79 

Vanadium IV 

3 f - *d 

3 f - 3 P 

7, 810. 79 

10, 611. 0 

(3d) 2 

9, 044. 00 

12. 930.4 

3 F - *G 

12, 031. 3 

18. 040. 0 

Chromium V 

3 f - *d 

3 f - 3 P 

9, 490. 51 

12. 647. 0 

(3d) 2 

10, 988.9 

15. 190. 0 

3 F - l G 

14, 618. 7 

21, 510. 0 


3 F - *D 

11, 170. 3 

14, 531. 

Manganese VI 
(3d) 2 

3 F - 3 P 

12. 933. 9 

17, 350. 7 

3 F - *G 

17, 206. 1 

24,706. 


3 F - *D 

12, 850.0 

16, 350. 

Iron VII 
(3d) 2 

3 f- 3 p 

14,878. 8 

19,455. 


3 F - *G 

19, 793. 5 

27, 790. 

Vanadium I 
(3d) 3 (4s) 2 

4 f- 4 p 

8, 363. 25 

9, 415. 24 
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Table 17- 

1 (con'd) 



Multiplet 

Separation 

Calculated 

Separation 

Observed 

Separation 

Chromium IV 
(3d) 3 

4 f- 4 p 

10. 308. 2 cm" 1 

13,816. cm" 

Manganese V 
(3d) 3 

4 f- 4 p 

12,253.2 

16,004. 

Iron VI 
(3d) 3 

4 f- 4 p 

14, 198. 1 

18,169. 

Cobalt VU 
(3d) 3 

4 f- 4 p 

16, 143. 1 

20. 271. 

Cobalt I 
(3d) 7 (4s) 2 

4 f- 4 p 

13,420. 1 

13. 068. 98 

Nickel I 
(3d) 8 (4s) 2 

3 f - 3 p 

14,684.4 

14, 604. 17 


The "observed separations" were obtained in all cases by averaging the fine structure 
splittings in the levels concerned. The experimental data were taken from Atomic 
Energy Levels, Charlotte E. Moore, Circular of the National Bureau of Standards 467, 

1949. - 

By definition, r(a) is the smaller, r(b) the larger of r and r'. The G k integral is ex¬ 
pressed as a similar sum. A table of L functions was made for n., n. and k values 
from zero to four, which is adequate to evaluate integrals involving electrons up to the 
3d shell. 

The ground state multiplet separations have been worked out for a number of 
elements. These results are given in Table 17-1. The multiplet separations for the 
(3d) 2 configurations reduce to linear functions of atomic number and are given by: 

3 F - 3 P * 5, 834. 8 Z 4 
3 F - l G = 7, 762. 1 Z. 

Z = 

3 1 4 " 

J F - D = 5, 039. 2 Z 4 

3 F - l S = 19,387. Z 4 
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There is good agreement between the slopes of the observed and theoretical 
separations as a function of Z. The (3d) 3 , *F - *P separation exhibits the same linear¬ 
ity in atomic number. The apparent accuracy in calculating the (3d) 7 and (3d)** multiplet 
separations is somewhat accidental. 
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18. AN ANALYTIC DETERMINATION OF ATOMIC ENERGY LEVELS IN MAGNESIUM 


Atomic energy levels for magnesium in various states of ionization and ex¬ 
citation were calculated using the determinantal method. ^ The experimental and cal¬ 
culated values for the configurations; Is 2 , ls 2 2s 2 2p*\ ls 2 2s 2 2p^3s 2 , ls 2 2s 2 2p^3s3p 
( 3 P and *P) are given in Table 18-1. For the closed shell configurations a single de- 


terminant was used. 

For the 3s3p energy level a 

secular equation 

corresponding to the 

spin-orbital degeneracy was jolved. 

Table 18-1 

Experimental 

Calculated 

Error 

Configuration 

value (cm* l ) 

value (cm - *) 

percentage 

is 2 

- 3.0011 x 10 7 

- 2. 9963 x 10 7 

. 15 

ls 2 2s 2 2p^ 

- 4.3759 x 10 7 

- 4. 3437 x 10 7 

. 74 

ls 2 2s 2 2p 6 3s3p 

3 P 

- 4.3931 x 10 7 

- 4. 3589 x 10 7 

. 78 

l P 

- 4.3944 x 10 7 

- 4. 3614 x 10 7 

. 78 


The one-electron wave functions used in the determinants were an orthonor¬ 
mal set formed from linear combinations of simple analytic shielding constant wave 
functions. ^ The normalized radial wave functions used are given in Table 18-2. 

Table 18-2 


1 0 80. 0403 e " 1L 7 r 

2 0 23.7035 e- 117 r - 36.7556 e' 3925 r 

2 1 35.2426 r e' 3> 925 r 

3 0 2. 57834 e' n ’ 7 r - 4. 25902 r e* 3 ' 925 r + . 354533 r 2 e'' 95 r 

3 1 3. 93672 r e' 3 ‘ 925 r -. 354539 r 2 e'• 95 r 
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19. THERMAL VIBRATIONS OF CU-ZN SYSTEM CRYSTALS 


The determination of the various branches of the normal vibration frequency 
versus reciprocal lattice vector contour map and the associated frequency spectra (num 
ber of normal frequency modes per unit frequency range) of the nuclear framework of 
crystals is of direct interest in the thermal diffuse scattering of x-rays, in restrahlen 
and Raman spectra interpretation, and in the calculation of such thermodynamic quan¬ 
tities as specific heats and coefficients of thermal expansion. There are three main 
aspects of the problem as treated on the Born-von Karman model: (l) the determina¬ 
tion of the generalized force constants between nuclei ("dynamical matrix" elements 
of Born^), which appear in the Taylor expansion of the nuclear potential energy (about 
the equilibrium cohesive energy value at which all first derivative or force terms are 
zero) as coefficients of the terms quadratic in the displacements of the nuclei from 
equilibrium; (2) the transformation of the coupled equations of motion of the nuclei into 
the normal coordinate form of linear, homogeneous, algebraic equations in the ampli¬ 
tudes of the normal coordinates; (3) the solution for the normal frequencies of the re¬ 
sulting consistency determinants of order three times the number of atoms per unit 
cell, and the calculation of the frequency spectrum either from this frequency plot or 
directly from these secular determinants. 

Inherent ir the method in the standard form, of course, is a neglect of quar- 
tic and higher order displacement terms in the potential energy, which would lead to 
anharmonic effects. ^ Step (2) above is equivalent to the well-known problem of diag- 
onalization of quadratic forms^ (here the simultaneous diagonalization of the kinetic 
energy quadratic terms in the velocities and the potential energy quadratic terms in the 
coordinates), and if the unit cell is simple this step presents little formal difficulty. 
Much attention has been focused on step (3), in which the sampling procedure of Black¬ 
man, ^ the moment method of Montroll, ^ and various graphical and mechanical 
methods for obtaining the frequency spectrum have all been used. The recent work of 
von Hove^ on the topology of the frequency contour map has provided a basic under¬ 
standing of the general features of the frequency spectrum. Quite often step (1) is not 
handled with as much care is the other two. Where real crystals are treated the force 
constants are usually determined from experimental values of the elastic constants. 

In order to reduce the number of independent parameters to the number of independent 
elastic constants and for general simplification, a number of assumptions are often 
made; namely, those of crystal isotropy; central forces; the potential being a sum of 
interaction energy terms between atoms taken only two at a time; and the vanishing of 
force constants for neighbors further away them second or third nearest (save for ionic 
crystals^). ^ 

In the present treatment of Cu-Zn system crystals particular attention is 
being paid to the problem of determining the force constants in as general and rigorous 
a way as possible. The symmetry properties of the crystals, as formulated in terms 
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(9) 

of the well-known theory of the representations of space groups, are of great help in 
simplifying this problem. ^ 10 ^ In obtaining the remaining independent force constants 
considerable use is to be made of the experimental results being obtained on the fre¬ 
quency contour maps of these crystals in the work of the x-ray group of Professor B. 

E. Warren at M. I.T. cn the temperature diffuse scattering of x-rays/ 1, ^ 
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